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ABSTRACT 

The thermal plasma of galaxy clusters lost most of its information on how structure formation 
proceeded as a result of dissipative processes. In contrast, non-equilibrium distributions of 
cosmic rays (CR) preserve the information about their injection and transport processes and 
provide thus a unique window of current and past structure formation processes. This infor- 
mation can be unveiled by observations of non-thermal radiative processes, including radio 
synchrotron, hard X-ray, and y-ray emission. To explore this, we use high-resolution simula- 
tions of a sample of galaxy clusters spanning a mass range of about two orders of magnitudes, 
and follow self-consistent CR physics on top of the radiative hydrodynamics. We model rela- 
tivistic electrons that are accelerated at cosmological structure formation shocks and those that 
are produced in hadronic interactions of CRs with ambient gas protons. We find that the CR 
proton pressure traces the time integrated non-equilibrium activities of clusters and is modu- 
lated by the recent dynamical activities. In contrast, the pressure of primary shock- accelerated 
CR electrons resembles current accretion and merging shock waves that break at the shallow 
cluster potential in the virial regions. The resulting synchrotron emission is predicted to be 
polarised and has an inhomogeneous and aspherical spatial distribution which matches the 
properties of observed radio relics. We propose a unified scheme for the generation of giant 
radio halos as well as radio mini-halos that naturally arises from our simulated synchrotron 
surface brightness maps and emission profiles. Giant radio halos are dominated in the centre 
by secondary synchrotron emission with a transition to the radio synchrotron radiation emit- 
ted from primary, shock- accelerated electrons in the cluster periphery. This model is able to 
explain the regular structure of radio halos by the dominant contribution of hadronically pro- 
duced electrons. At the same time, it is able to account for the observed correlation of mergers 
with radio halos, the larger peripheral variation of the spectral index, and the large scatter 
in the scaling relation between cluster mass and synchrotron emission. Future low-frequency 
radio telescopes (LOFAR, GMRT, MWA, LWA) are expected to probe the accretion shock re- 
gions of clusters and the warm-hot intergalactic medium, depending on the adopted model for 
the magnetic fields. The hadronic origin of radio halos can be scrutinised by the detection of 
pion-decay induced y-rays following hadronic CR interactions. The high-energy y-ray emis- 
sion depends only weakly on whether radiative or non-radiative gas physics is simulated due 
to the self-regulated nature of the CR cooling processes. Our models predict a y-ray emission 
level that should be observable with the GLAST satellite. 

Key words: cosmology: large-scale structure of Universe, galaxies: cluster: general, mag- 
netic fields, cosmic rays, radiation mechanisms: non-thermal 



1 INTRODUCTION 
1.1 Motivation 

A substantial number of luminous X-ray clusters show diffuse 
large-scale radio emission. Generally these radio phenomena can 



* e-mail: pfrommer@cita.utoronto.ca (CP); ensslin@mpa- 
garching.mpg.de (TAE); voLk:er@mpa-garching.mpg.de (VS) 



be divided into two categories that differ morphologically, in their 
degree of polarisation, as well as in their characteristic emission 
regions with respect to the cluster halo. The large-scale "radio 
relic" or "radio gischt" emission (Kempner et al. 2004), that has 
a high degree of polarisation, is irregularly shaped and occurs at 
peripheral cluster regions, can be attributed to merging or accretion 
shock waves as proposed by Ensslin et al. (1998). Prominent ex- 
amples for large scale "radio relic" emission have been observed 



© 2003 RAS 



2 C. Pfrommer, T. A. Enfilin, V, Springel 



in Abell 3667 (Rottgering et al. 1997), Abell 3376 (Bagchi et al. 
2006), and Abell 2256 (Bridle & Fomalont 1976; Masson & Mayer 
1978; Bridle et al. 1979; Rottgering et al. 1994; Clarke & EnBlin 
2006). In contrast, the origin of "cluster radio halos" that resemble 
the regular morphology of the X-ray emitting intra-cluster plasma 
is not understood to date. Prominent examples for "radio halo" 
emission can be obtained from Giovannini et al. (1999) and include 
the Coma cluster (Kim et al. 1989; Deiss et al. 1997) and the galaxy 
cluster IE 0657-56 (Liang et al. 2000). In principle, observations 
of non-thermal cluster phenomena could provide an independent 
and complementary way of studying the growth of structure in our 
Universe and could shed light on the existence and the properties of 
the warm-hot intergalactic medium (WHIM), provided the underly- 
ing processes are understood. Sheets and filaments are predicted to 
host this WHIM with temperatures in the range 10^ K < T < 10^ K 
whose evolution is primarily driven by shock heating from gravita- 
tional perturbations breaking on mildly nonlinear, non-equilibrium 
structures (Hellsten et al. 1998; Cen & Ostriker 1999; Dave et al. 
2001; Furlanetto & Loeb 2004; Kang et al. 2005). Once a clus- 
ter has virialised, the thermal plasma lost most information on 
how the formation proceeded due to the dissipative processes driv- 
ing the plasma towards a Maxwell-Boltzmann momentum distribu- 
tion that is characterised by its temperature only. In contrast, non- 
equilibrium distributions of cosmic rays preserve the information 
about their injection and transport processes much better, and thus 
provide a unique window of current and past structure formation 
processes. 

The information about these non-equilibrium processes is en- 
coded in the spectral and spatial distribution of cosmic ray elec- 
trons and protons. Radiative loss processes of these non-thermal 
particle distributions produce characteristic radio synchrotron, hard 
X-ray inverse Compton, and hadronically induced y-ray emis- 
sion. Suitably combining various non-thermal emission processes 
will allow us to infer the underlying non-equilibrium processes of 
clusters as well as to improve our knowledge about fundamental 
plasma physics. The upcoming generation of radio, hard X-ray, 
and y-ray instruments opens up the extragalactic sky in unexplored 
wavelength ranges: low-frequency radio arrays (GMRT^, LOFAR^, 
MWA^, LWA^), the future hard X-ray satellite missions NuSTAR^ 
and Simbol-X, and y-ray instruments (the GLAST^ satellite and 
imaging atmospheric Cerenkov telescopes H.E.S.SJ , MAGIC^, 
VERITAS^, and CANGAROO^^) will allow us to probe non-thermal 
cluster physics with a multi-faceted approach. There have been pi- 
oneering eff'orts to simulate the non-thermal emission from clus- 
ters by numerically modelling discretised cosmic ray (CR) energy 
spectra on top of Eulerian grid-based cosmological simulations 
(Miniati 2001; Miniati et al. 2001a,b; Miniati 2002, 2003). How- 
ever, these approaches neglected the hydrodynamic pressure of the 
CR component, were unable to resolve the observationally accessi- 
ble, dense central regions of clusters, and they neglected dissipative 
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Figure 1. Overview over the relevant physical processes in galaxy clus- 
ters. The right side shows the interplay of different physical processes high- 
lighting the interplay of the energy reservoirs of the thermal plasma and 
cosmic ray protons (shown in blue) while the left side shows observables 
that inform about the properties of clusters and their dynamical state. Gain 
processes are denoted in green, while loss or redistribution processes are 
denoted in red. 

gas physics including radiative cooling, star formation, and super- 
nova feedback. To allow studies of the dynamical eff'ects of CR 
protons in radiatively cooling galactic and cluster environments, we 
have developed a CR proton formalism that is based on smoothed 
particle hydrodynamical representation of the equations of motion. 
The emphasis is given to the dynamical impact of CR protons on 
hydrodynamics, while simultaneously allowing for the important 
CR proton injection and loss processes in a cosmological setting 
(EnBlin et al. 2007; Jubelgas et al. 2007; Pfrommer et al. 2006). 
This enables us to account for the pressure support provided by CR 
protons to the plasma of clusters of galaxies. A substantial CR pro- 
ton pressure contribution might have a major impact on the prop- 
erties of the intra-cluster medium (ICM) and could modify thermal 
cluster observables such as the X-ray emission and the Sunyaev- 
Zel'dovich (SZ) eff^ect (Pfrommer et al. 2007). 

1.2 Cosmic ray physics 

We give a short and simplified overview over the relevant physical 
processes in galaxy clusters in Fig. 1 before introducing diff'erent 
CR populations that are relevant for the non-thermal emission from 
clusters. The upper central part of Fig. 1 shows standard processes 
which are usually considered in simulations. Radiative cooling of 
the gas leads eventually to star formation in the densest regions 
that exceed a certain density threshold. This happens in the cen- 
tral cluster regions and within substructures leading to individual 
galaxies. Once the nuclear energy has been used up, massive stars 
explode in supernovae that drive strong shock waves into the am- 
bient interstellar medium (ISM) which resupply thermal and tur- 
bulent energy. On larger scales, structure formation shock waves 
dissipate gravitational energy associated with hierarchical cluster- 
ing into thermal energy of the gas, thus supplying the ICM with 
entropy and thermal pressure support. There are three main observ- 
ables associated with these processes: the hot ICM emits thermal 
bremsstrahlung radiation with an X-ray luminosity that depends on 
the square of the electron number density. The amplitude of the 
Sunyaev-Zel'dovich eff'ect (Sunyaev & Zeldovich 1972) depends 
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Figure 2. Schematic overview over non-thermal radiative processes in 
galaxy clusters. Various gravitational and non-gravitational energy sources 
(shown in red) are able to accelerate relativistic particle populations (shown 
in blue) by means of different plasma processes (shown in green). Non- 
thermal cluster observables (shown in yellow) are tracers of these cosmic 
ray populations: any cosmic ray electron population can emit radio syn- 
chrotron radiation as well as inverse Compton emission that extends from 
the X-ray into the y-ray regime. In contrast, the characteristic spectral signa- 
ture accompanying y-ray emission from hadronic cosmic ray interactions is 
a unique sign of a cosmic ray proton population in the intra-cluster plasma. 



on the pressure of the thermal electron population integrated along 
the line-of- sight through the cluster. Finally, galaxy spectra probe 
directly the stellar populations of intra-cluster galaxies and indi- 
rectly the cluster's potential through their velocity dispersion (for 
reviews see Sarazin 1988; Voit 2005). 

The lower part of Fig. 1 sketches the cosmic ray physics 
within clusters. CR protons behave differently compared to the 
thermal gas. Their equation of state is softer, they are able to travel 
actively over macroscopic distances, and their energy loss time- 
scales are typically larger than the thermal ones. Besides ther- 
malization, collisionless shocks are also able to accelerate ions 
of the high-energy tail of the Maxwellian through diffusive shock 
acceleration (for reviews see Drury 1983b; Blandford & Eichler 
1987; Malkov & O'C Drury 2001). These energetic ions are re- 
flected at magnetic irregularities through magnetic resonances be- 
tween the gyro-motion and waves in the magnetised plasma and 
are able to gain energy in moving back and forth through the 
shock front. This acceleration process typically yields a CR proton 
population with a power-law distribution of the particle momenta. 
CRs are accelerated on galactic scales through supernova shocks 
while they are injected by structure formation shock waves on even 
larger scales up to tens of Mpc. So far, we have neglected feed- 
back from active galactic nuclei (AGN) in our simulations despite 
its importance (for first numerical simulations of thermal 'radio- 
mode' feedback within cosmological simulations, see Sijacki & 
Springel 2006). Gravitational energy associated with the accretion 
onto super-massive black holes is converted into large-scale jets 
and eventually dissipated into thermal and CR energy. 

The ubiquitous cosmic magnetic fields couple the otherwise 
dynamically independent ingredients like the ICM plasma, and the 
CR gas into a single, however complex fluid. Magnetic fields pre- 
vent charged relativistic particles to travel macroscopic distances 
with their intrinsic velocity close to the speed of light. Instead, 
the particles gyrate around, and travel slowly along magnetic field 



lines. Occasionally, they get scattered on magnetic irregularities. 
On macroscopic scales, the transport can often be described as a 
diffusion process that redistributes the CR energy density macro- 
scopically provided the gyro-radius of charged relativistic particles 
can be regarded to be small. Thus, the diffusive CR transport in 
tangled magnetic fields effectively confines the CRs with energies 
£^ < 2 X 10^ GeV to clusters and yields to CR proton lifetimes of 
the order of the Hubble time (Volk et al. 1996; Ensslin et al. 1997; 
Berezinsky et al. 1997; Colafrancesco & Blasi 1998), long enough 
to diffuse away from the production site and to maintain a space- 
filling distribution over the cluster volume. Thermal heat conduc- 
tion is an analogous process that reallocates the thermal energy of 
the ICM. 

The CR energy reservoir suffers two main loss processes: (1) 
CR energy is transferred into the thermal energy reservoir through 
individual electron scatterings in the Coulomb field of the CR par- 
ticle as well as by small momentum transfers through excitations of 
quantised plasma oscillations. We refer to the sum of both effects 
as Coulomb losses (Gould 1972a). (2) Provided the CR momentum 
exceeds the threshold /? ^ 0.8 GeV/c for the hadronic reaction with 
ambient protons, they produce pions which decay into secondary 
electrons, positrons, neutrinos, and y-ray s: 

n- fi- + y,u/yiu ^ e- + Ve/ve + yiu + yiu 

Only CR protons above this kinematic threshold are therefore vis- 
ible through their decay products via radiative processes, making 
them directly observationally detectable. As shown in Fig. 2, these 
secondary relativistic electrons and positrons can emit a halo of 
radio synchrotron emission in the presence of ubiquitous intra- 
cluster magnetic fields (Dennison 1980; Vestrand 1982; Blasi & 
Colafrancesco 1999; Dolag & EnBlin 2000; Miniati et al. 2001b; 
Pfrommer & EnBlin 2003, 2004a,b; Marchegiani et al. 2007) as 
well as inverse Compton emission by scattering photons from the 
cosmic microwave background into the hard X-ray and y-regime.^^ 
Future y-ray satellites should be able to detect the associated 
hadronically induced y-ray emission resulting from neutral pion 
decay and allow unambiguous conclusions on the parent CR popu- 
lation in clusters. 

Structure formation shocks can also directly accelerate so- 
called primary CR electrons giving rise to an irregularly shaped 
radio and inverse Compton morphology due to the comparatively 
short synchrotron lifetimes of CR electrons of r ^ 10^ yr. To com- 
plicate this picture even more, there are other processes that ac- 
celerate relativistic electrons. Re- acceleration processes of 'mildly' 
relativistic electrons (y ^ 100 - 300) that are being injected over 
cosmological timescales into the ICM by sources like radio galax- 
ies, supernova remnants, merger shocks, or galactic winds can pro- 
vide an efficient supply of highly-energetic CR electrons. Owing 
to their long lifetimes of a few times 10^ years these 'mildly' rel- 
ativistic electrons can accumulate within the ICM (Sarazin 2002), 
until they experience continuous in- situ acceleration either via in- 
teractions with magneto-hydrodynamic waves, or through turbulent 
spectra (Jaffe 1977; Schlickeiser et al. 1987; Brunetti et al. 2001; 
Ohno et al. 2002; Brunetti et al. 2004; Brunetti & Lazarian 2007). 
This gives rise to a third population of re-accelerated CR electrons 
that also contributes to the observed radio and inverse Compton 
emission. Since the distribution of magnetic field strengths with 



In the following, we use the term secondary CR electrons synonymously 
for the likewise produced electrons and positrons. 
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cluster radius is also not well known, radio synchrotron emission 
alone has limited predictive power. Unfortunately, the conceptually 
simpler inverse Compton emission is hard to observe because of 
the strong radiation background in the soft and hard X-ray regime. 

Nevertheless, there seems to be growing evidence for an ex- 
cess of hard X-ray emission compared to the expected thermal 
bremsstrahlung in a number of clusters that is based on obser- 
vations with instruments on board five diff'erent X-ray satellites. 
Prominent examples include the Coma cluster (Rephaeli et al. 
1999; Fusco-Femiano et al. 1999; Rephaeli & Gruber 2002; Fusco- 
Femiano et al. 2004, 2007;^^ Eckert et al. 2007; using the Rossi 
X-ray Timing Explorer (RXTE), BeppoSAX, and INTEGRAL) , 
the "Bullet" cluster lES 0657-558 (Petrosian et al. 2006, using 
RXTE), Abell 2256 (Rephaeli & Gruber 2003, using RXTE), the 
Ophiuchus cluster (Eckert et al. 2007, using INTEGRAL), and the 
Perseus cluster (Sanders et al. 2005; Molendi 2007, using Chandra 
and XMM-Newton). The currently favoured theoretical explanation 
of inverse Compton radiation by CR electrons faces problems since 
the magnetic field estimates inferred by combining synchrotron and 
inverse Compton emission are much smaller than those derived 
from Faraday rotation measurements (cf. Pfrommer 2007, for an 
extensive discussion). It has been proposed that a fraction of the 
diff'use cosmological y-ray background radiation originates from 
the same processes (Loeb & Waxman 2000; Miniati 2002; Reimer 
et al. 2003; Berrington & Dermer 2003; Kuo et al. 2005). 

This paper studies directly the CR related multi-frequency ob- 
servables aiming at understanding the cluster radio halo emission. 
In a companion paper, we study the interplay of thermal gas and 
CRs and their eff'ect on the observables of the thermal gas such as 
X-ray emission and the Sunyaev-Zel'dovich eff'ect (Pfrommer et al. 
2007, hereafter Paper I). For consistency reasons with that paper, 
we scale cluster masses and length units with the dimensionless 
Hubble constant, h, where Hq = lOO/zkms"^ Mpc~^ However, 
non-thermal surface brightness and luminosities (for radio, hard 
X-ray, and y-ray emission) are reported in units of the currently 
favoured Hubble constant, hjo, where Hq = lOhjo km s~^ Mpc~^ 
since primary and secondary emission processes have a diff'erent 
scaling with h. We derive cluster scaling relations for non-thermal 
observables and compare our results to observations in our follow- 
up paper (Pfrommer 2007, hereafter Paper III). 

The outline of the paper is as follows. Section 2 describes 
our methodology including the general setup of the simulations, 
our cluster sample, the diff'erent physical processes which we sim- 
ulated, and highlights important properties of radiative processes 
considered in this work. In Sect. 3, we present and interpret the 
results on the line-of-sight projections and emission profiles of 
the diff'erent non-thermal emission mechanisms, and correlations 
of various non-thermal emission processes with the thermal X-ray 
emission. We compare our results to previous findings in the litera- 
ture and point out future theoretical work that is needed to comple- 
ment this work (cf. Sect. 4). The conclusions are drawn in Sect. 5. 
Appendix A describes the modelling of the primary and secondary 
CR electron population while Appendix B describes the formulae 
of the non-thermal emission mechanisms ranging from radio syn- 
chrotron radiation, inverse Compton emission, as well as hadroni- 
cally induced y-ray emission. 



The results of these papers have been challenged by an analysis that 
takes into account all systematic uncertainties in the critical parameters in- 
cluding the choice of a source-free background field and the modelling of 
the thermal model for the ICM (Rossetti & Molendi 2004, 2007). 



2 METHODOLOGY 

2.1 General approach 

We follow the CR proton pressure dynamically in our simulations 
while taking into account all relevant CR injection and loss terms 
in the ICM, except for a possible production of CR protons by 
AGN. In contrast, we model the CR electron population in a post- 
processing step since it does not modify the hydrodynamics owing 
to its negligible pressure contribution. In this paper, we concen- 
trate on three observationally motivated wave-bands. (1) Radio syn- 
chrotron emission between 15 MHz and 1.4 GHz, (2) non- thermal 
hard X-ray emission at energies Ey > 10 keV, and (3) y-ray emis- 
sion at energies E^ > 100 MeV. Studying our simulated radio syn- 
chrotron maps and emission profiles, we develop a unified scheme 
for the generation of the diffiise large scale radio emission of clus- 
ters such as giant radio halos, mini-halos, as well as the radio relic 
emission. 

2.2 The simulations 

2. 2. 1 General setup and cluster sample 

This section provides a short overview of the simulations and phys- 
ical models used. Further details can be found in Paper I. All simu- 
lations were performed using the "concordance" cosmological cold 
dark matter model with a cosmological constant (ACDM). The cos- 
mological parameters of our model are: Qm = ^dm + = 0.3, 
= 0.039, = 0.7, h = 0.7, w = 1, and erg = 0.9. Here, 
denotes the total matter density in units of the critical density 
for geometrical closure today, pcnt = 3Hl/{SnG). Q.^ and Q.a de- 
note the densities of baryons and the cosmological constant at the 
present day. The Hubble constant at the present day is parametrised 
as Hq = 100 /z km s~^Mpc~^ while n denotes the spectral index of 
the primordial power-spectrum, and erg is the rms linear mass fluc- 
tuation within a sphere of radius 8 h~^Mpc extrapolated to z = 0. 

Our simulations were carried out with an updated and ex- 
tended version of the distributed-memory parallel TreeSPH code 
GADGET-2 (Springel 2005; Springel et al. 2001) that includes self- 
consistent cosmic ray physics (EnBlin et al. 2007; Jubelgas et al. 
2007; Pfrommer et al. 2006). Gravitational forces were computed 
using a combination of particle-mesh and tree algorithms. Hydro- 
dynamic forces are computed with a variant of the smoothed parti- 
cle hydrodynamics (SPH) algorithm that conserves energy and en- 
tropy where appropriate, i.e. outside of shocked regions (Springel 
& Hernquist 2002). 

We have performed high-resolution hydrodynamic simula- 
tions of the formation of 14 galaxy clusters. The clusters span a 
mass range from 5x10^^ h~^ M© to 2 x 10^^ h~^ M© and show a va- 
riety of dynamical states ranging from relaxed cool core clusters to 
violent merging clusters (cf. Table 1). The clusters have originally 
been selected from a low-resolution dark-matter-only simulation 
(Yoshida et al. 2001). Using the 'zoomed initial conditions' tech- 
nique (Katz & White 1993), the clusters have been re-simulated 
with higher mass and force resolution by adding short-wavelength 
modes within the Lagrangian regions in the initial conditions that 

Our 'unified scheme' unifies apparently difi'erent difi'use radio phenom- 
ena in clusters (giant relics, halos, and mini-halos) with a simple and physi- 
cally motivated model. The 'unified scheme' should not be confused with a 
'complete model' and we want to point out that we have not accounted for 
all possible CR processes that could be of interest in the context of cluster 
physics. 
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Table 1. : Cluster sample 



Cluster 


sim.'s 


dyn. state*^^^ 


, ,(2) 
^200 

\h-^ Mo] 


^200 

\h~^ Mpc] 


'^^ 200 
[keV] 


1 


28a 


CC 


1.8 X 10^^ 


2.0 


13.1 


2 




cc 


1.3 X 10^^ 


1.8 


10.6 


3 


s72a 

t> ' 


PostM 


1.1 X 10^^ 


1.7 


9.4 


4 


251 


CC 


1.1 X 10^^ 


1.7 


9.4 


5 




M 


3.7 X 10^"^ 


1.2 


4.7 


6 


g72b 


M 


1.5 X 10^4 


0.87 


2.4 


7 


glc 


M 


1.4 X 10^4 


0.84 


2.3 


8 


g8b 


M 


l.Ox IQi^^ 


0.76 


1.9 


9 


gld 


M 


9.2 X 10^3 


0.73 


1.7 


10 


g676 


CC 


8.8 X 10^3 


0.72 


1.7 


11 


g914 


CC 


8.5 X 10^3 


0.71 


1.6 


12 


gle 


M 


6.4 X 10^3 


0.65 


1.3 


13 


g8c 


M 


5.9 X 10^3 


0.63 


1.3 


14 


g8d 


PreM 


5.4 X 10^3 


0.61 


1.2 


Notes: 



(1) The dynamical state has been classified through a combined criterion 
invoking a merger tree study and the visual inspection of the X-ray bright- 
ness maps. The labels for the clusters are M-merger, PostM-post merger 
(slightly elongated X-ray contours, weak cool core region developing), 
PreM-pre-merger (sub-cluster already within the virial radius), CC-cool 
core cluster with extended coohng region (smooth X-ray profile). 

(2) The virial mass and radius are related by Ma(z) = ^tt Apcrit(^)^^5 
where A = 200 denotes a multiple of the critical overdensity Pcrit(^) = 
?>H{zfl{^nG). 

(3) The virial temperature is defined by kT^ = GM/^/dmp/(2R/^), where ji 
denotes the mean molecular weight. 



will evolve later-on into the structures of interest (Paper I). We re- 
simulated three isolated clusters (cluster 4, 10, and 11) and three 
super-cluster regions which are each dominated by a large cluster 
(cluster 1, 2, and 3) and surrounded by satellite clusters (cluster 5 
to 9 and 12 to 14). In high-resolution regions, the dark matter par- 
ticles had masses of mbm = 1.13 x 10^ M© and SPH particles 
had mgas = 1.7 X 10^ M© so each individual cluster is resolved 
by 8 X 10"^ to 4 X 10^ particles, depending on its final mass. The 
SPH densities were computed from 48 neighbours, allowing the 
SPH smoothing length to drop at most to half of the value of the 
gravitational softening length of the gas particles. This choice of 
the SPH smoothing length leads to our minimum gas resolution 
of approximately 8 x 10^ M©. For the initial redshift we chose 
1 + ^init = 60. The gravitational force softening was of a spline form 
(e.g., Hernquist & Katz 1989) with a Plummer equivalent softening 
length that is assumed to have a constant comoving scale down to 
z = 5, and a constant value of 5 h~^kpc in physical units at later 
epochs. 

We analysed the clusters with a halo-finder based on spherical 
overdensity followed by a merger tree analysis in order to get the 
mass accretion history of the main progenitor. We also produced 
projections of the X-ray emissivity at redshift z = in order to 
get a visual impression of the cluster morphology. The dynamical 
state of a cluster is defined by a combined criterion: (i) if the clus- 
ter did not experience a major merger with a progenitor mass ratio 
1:3 or larger after z = 0.8 (corresponding to a look-back time of 
~ 5 Gyr) and (ii) if the visual impression of the cluster's X-ray 
morphology is relaxed, it was defined to be a cool core cluster. The 
spherical overdensity definition of the virial mass of the cluster is 
given by the material lying within a sphere centred on a local den- 
sity maximum, whose radial extend Ra is defined by the enclosed 



Table 2. : Different physical processes in our simulations: 



Simulated physics^^^ simulation models^^-*: 
Si S2 S3 



thermal shock heating 








radiative cooling 




/ 


/ 


star formation 




/ 


/ 


Coulomb CR losses 


>/ 


/ 


/ 


hadronic CR losses 




/ 




shock-CRs 






/ 


supemova-CRs 









Notes: 

(1) This table serves as an overview over our simulated models. The first 
column shows the simulated physics and the following three columns show 
our diff'erent simulation models with varying gas and cosmic ray physics. 
Model SI models the thermal gas non-radiatively and includes CR physics, 
while the models S2 and S3 use radiative gas physics with diff"erent variants 
of CR physics. 



threshold density condition M(< /?A)/(47ri?^/3) = Pthres- We chose 
the threshold density Pthres(^) = Apcrit(z) to be a multiple A = 200 
of the critical density of the universe Pcritfe) = 3//(z)^/(8;7rG). We 
assume a constant A = 200 although some treatments employ a 
time- varying A in cosmologies with Q.^ ^ 1 (Eke et al. 1996). In 
the reminder of the paper, we use the terminology /?vir instead of 

-^200- 

2.2.2 The models 

For each galaxy cluster we ran three diflferent simulations with 
varying gas and cosmic ray physics (cf. Table 2). The first set of 
simulations used non-radiative gas physics only, i.e. the gas is trans- 
ported adiabatically unless it experiences structure formation shock 
waves that supply the gas with entropy and thermal pressure sup- 
port. Additionally we follow cosmic ray (CR) physics including 
adiabatic CR transport processes, injection by cosmological struc- 
ture formation shocks with a Mach number dependent acceleration 
scheme, as well as CR thermalization by Coulomb interaction and 
catastrophic losses by hadronic interactions (model SI). The sec- 
ond set of simulations follows the radiative cooling of the gas, star 
formation, supernova feedback, and a photo-ionising background 
(details can be found in Paper I). As before in model SI, we ac- 
count for CR acceleration at structure formation shocks and allow 
for all CR loss processes (model S2). The last set of simulations ad- 
ditionally assumes that a constant fraction ^sn = ^CR,inj /f^diss = 0.3 
of the kinetic energy of a supernova ends up in the CR population 
(model S3), which is motivated by TeV y-ray observations of a su- 
pernova remnant that find an energy fraction of ^sn - 0. 1 -0.3 when 
extrapolating the CR distribution function (Aharonian et al. 2006). 
We choose a maximum value for this supernova energy efficiency 
owing to the large uncertainties and our aim to bracket the realis- 
tic case with the two radiative CR simulations. Generally, we use 
model S2 as our standard case which is a conservative choice for 
the CR pressure and explore how the physics of the other models 
change the resulting non-thermal cluster observables. In this work, 
we don't account for AGN sources of cosmic rays, but that this will 
be addressed in upcoming work (Sijacki et al. 2007). 

Radiative cooling was computed assuming an optically thin 
gas of primordial composition (mass-fraction of Xh = 0.76 for hy- 
drogen and 1 - Xh = 0.24 for helium) in coUisional ionisation 
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equilibrium, following Katz et al. (1996). We have also included 
heating by a photo-ionising, time-dependent, uniform ultraviolet 
(UV) background expected from a population of quasars (Haardt 
& Madau 1996), which reionises the Universe at z - 6. Star forma- 
tion is treated using the hybrid multiphase model for the interstellar 
medium introduced by Springel & Hernquist (2003). In short, the 
ISM is pictured as a two-phase fluid consisting of cold clouds that 
are embedded at pressure equilibrium in an ambient hot medium. 
The clouds form from the cooling of high density gas, and rep- 
resent the reservoir of baryons available for star formation. When 
stars form, the energy released by supernovae heats the ambient hot 
phase of the ISM, and in addition, clouds in supernova remnants 
are evaporated. These effects establish a tightly self-regulated sub- 
resolution model for star formation in the ISM. 

Cosmic ray physics was computed by using a new formula- 
tion that follows the most important injection and loss processes 
self-consistently while accounting for the CR pressure in the equa- 
tions of motion (EnBlin et al. 2007; Jubelgas et al. 2007; Pfrom- 
mer et al. 2006). We refer to these papers for a detailed descrip- 
tion of the formalism, providing here only a short summary of the 
model. In our methodology, the non- thermal cosmic ray population 
of each gaseous fluid element is approximated by a simple power 
law spectrum in particle momentum, characterised by an ampli- 
tude, a low-momentum cut-off, and a fixed slope a = 2.3. This 
choice is justified by taking the mean of the Mach number distribu- 
tion weighted by the dissipated energy at shocks which is closely 
related to the spectral index of a CR power-law distribution (Pa- 
per I). Adiabatic CR transport processes such as compression and 
rarefaction, and a number of physical source and sink terms which 
modify the cosmic ray pressure of each particle are modelled. The 
most important sources considered are injection by supernovae (in 
our radiative simulations) and diffusive shock acceleration at cos- 
mological structure formation shocks, while the primary sinks are 
thermalization by Coulomb interactions, and catastrophic losses by 
hadronic interactions. 

2.3 Essentials of radiative processes 

We are interested in the non-thermal emission of the inter-galactic 
medium at radio frequencies (y > 10 MHz) as well as at hard X- 
ray/y-ray energies (Ey > 10 keV). This non- thermal emission is 
generated by CR electrons with energies > GeV as can be read- 
ily inferred from the classical synchrotron and inverse Compton 
formulae, 

3eB . B ( y \^ 

= 2;^^^^^«^Hg(i^)' 

hvic = ^/iyinit/-90keV^ (-^f, (2) 

where e denotes the elementary charge, h the Planck constant, c 
the speed of light, mg the electron mass, the particle kinetic energy 
EKtUqC^) = 7 - 1 is defined in terms of the Lorentz factor y, and 
B = ^{B^) is the rms of the magnetic vector field B. We chose 
CMB photons /zvcmb - 0-66 meV as source for the inverse Comp- 
ton emission using Wien's displacement law. Thus, the same CR 
electron population seen in the radio band via synchrotron emission 
can be observed in the hard X-ray regime through the IC process. 

2. 3. 1 Synchrotron and IC emission from equilibrium spectra 

The synchrotron and IC emissivities of an equilibrium distribu- 
tion of CR electrons is characterised by two distinctive regimes. 
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Figure 3. The synchrotron and inverse Compton (IC) emissivity of an equi- 
librium distribution of CR electrons is shown for various spectral indices ay. 
The normalisation is given by the respective emissivities at the equivalent 
magnetic field strength of the energy density of the CMB, 5cmb = 3.24 (oG. 
The weak field regime on the left-hand side is characterised by the domi- 
nant inverse Compton (IC) emission while the strong field regime on the 
right side has the synchrotron emission as the dominant electron cooling 
channel. 

The weak field regime is characterised by the dominating IC emis- 
sion while the strong field regime has the synchrotron emission as 
the dominant electron cooling channel. Using the formulae for the 
equilibrium distribution of CR electrons (Eqns. (A35) or (A48)), 
the synchrotron and IC emissivities depend on the magnetic field 
strength as follows: 

jv ^ and 7ic oc , (3) 

+ ^CMB + ^CMB 

where ay denotes the synchrotron spectral index that is defined by 
jy oc y~^^, and Sb - B^/i^n). Figure 3 shows these emissivities, 
normalised by respective emissivities at the equivalent magnetic 
field strength of the CMB energy density at z = 0, Bqmb = 3.24 |xG. 
In the IC dominated electron cooling regime leftwards of Bqmb^ the 
synchrotron emissivity quickly decreases relative to the IC emis- 
sivity, showing thus a strong dependence on the magnetic field 
strength. There is an interesting twist associated with the differ- 
ent spectral indices of the synchrotron emission in clusters. Paper I 
shows that the characteristic shock strength increases as one moves 
outwards from the cluster centre due to the decrease of the sound 
velocity in combination with the shallower peripheral cluster po- 
tential. CR acceleration crucially depends on the shock strength 
according to Eqn. (A 19) predicting steep CR spectra at the cen- 
tre that flatten on average towards the cluster periphery. Since the 
magnetic field has a decreasing profile with radius (Eqn. 5), the 
synchrotron emission of clusters should qualitatively be given by 
the upper envelope of the family of emissivity curves labelled by 
different spectral indices of Fig. 3. This simplified picture assumes 
that the electron spectra are dominated by injection and are nei- 
ther qualitatively modified by CR transport processes such as CR 
diffusion nor by the hadronically injected electron spectra. 

2.3.2 Comparison of inverse Compton and 'j-ray spectra 

Figure 4 compares the spectral distribution of the pion-decay in- 
duced y-ray emission (broken power-laws) with the inverse Comp- 
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Ey [GeV] 

Figure 4. Spectral distribution of the differential y-ray flux dTy/dEy (up- 
per panel) and the integrated y-ray flux Ty(Ey, oo) (bottom panel) for dif- 
ferent spectral indices. Shown are inverse Compton spectra from a sec- 
ondary CR electron population and pion decay induced y-ray spectra (bro- 
ken power-laws) both resulting from hadronic CR proton interactions. The 
model calculations assume a normalisation for the CR proton distribution 
of Cp = 10~^^ cm~^, a nucleon density of = 10~^ cm~^ (for Ay), 
and put the fiducial cluster with mass Mduster = 10^^ Mq/H at a distance 
D = 100Mpc//z to derive Ty. 



ton emission both resulting from hadronic CR proton interactions. 
Note, that the relative normalisation of both emission components 
is governed by hadronic physics and does not depend on the gas 
and CR proton number densities. For our choice of the CR proton 
spectral index of ofp = 2.3, the ratio of pion-decay to secondary IC 
emission in the energy range Ey > 100 MeV can be readily inferred 
to be Ty/Tic = 20. The asymptotic behaviour for the y-ray number 
flux of both emission components at high energies is given by 

Ty oc ^7^^' and Tic oc = ^^p^'. (4) 

Assuming a spectral CR index of ofp = 2 yields the same asymptotic 
behaviour while increasing results in a shallower decline for Tic 
with energy such that eventually the IC component will surpass the 
pion decay emission. This however is well above the energy range 
Ey \0 TeV that is of interest for imaging air Cerenkov tele- 
scopes. For an energy range Ey < 1 MeV the secondary IC emis- 
sion dominates the hadronically induced channel. In contrast to the 
secondary IC emission, the IC emission level of primary CR elec- 
trons depends on the dynamical activity of the region, in particular 
on the shock strength and the associated amount of dissipated en- 
ergy. Comparing the primary to the secondary IC emission will be 
one goal of this paper. 



3 RESULTS 

3.1 Cluster environment and cosmic ray pressure 

Figure 5 shows the region around the largest cluster with merg- 
ing activity in our sample (M ^ lO^^h'^ M©, with the identifier 
g72a) in our simulation with radiative gas physics, star formation, 
and CRs from structure formation shocks only (model S2). This 
galaxy cluster experienced a large merger with a mass ratio of 
'Wmerger/'Wprogenitor = 0.3 at rcdshift z = 0.25 prcccdcd by a minor 
merger mass ratio of mmerger/'Wprogenitor = 0.1 at z = 0.4. The en- 
vironment is dominated by the large central post-merging cluster 
and surrounded by smaller satellite clusters and groups. The line-of 
sight average of the density shows the suitably normalised quantity 
1 + ^gas = p/(^h Pcrit). While the ICM of the central massive clus- 
ter reaches central temperatures above that of the virial temperature 
of kTyij. = 9.4 keV, the surrounding warm-hot intergalactic medium 
(WHIM) acquires temperatures of kT ~ (10~^-10~^) keV. The spa- 
tial distribution of shock strengths can be studied best by looking at 
the Mach numbers weighted by the energy dissipation rate at struc- 
ture formation shocks (represented by the colour hue in the bottom 
left panel of Fig. 5). The brightness scales logarithmically with the 
projected dissipation rate <6:diss>ios- Within this super-cluster region 
most of the energy is dissipated in weak internal shocks with Mach 
numbers At < 2 which are predominantly central flow shocks or 
merger shock waves traversing the cluster centre. Collapsed cosmo- 
logical structures are surrounded by several quasi- spherical shells 
of external shocks with successively higher Mach numbers, but 
they play only a minor role in the energy balance of thermaliza- 
tion as can be inferred by its dim brightness. Clearly visible are 
spherical shells of shocks at diflferent radii from the cluster centre. 
Two distinct outgoing shock waves at distances of 2 and 3 Mpc 
to the cluster centre are visible that are triggered by the merger, 
followed by shells of stronger shocks further outwards. This pic- 
ture is dramatically changed if we perform the weighting of the 
Mach numbers with the energy injection rate of CR protons, ^cr 
(shown in colour while the brightness displays the logarithm of the 
CR proton energy injection rate, bottom right side). Only CR pro- 
tons with a dimensionless momentum q = /3y > 30 (E > 30 GeV) 
have been considered for calculating the CR energy density since 
lower energetic CR protons are not detectable at radio frequencies 
y > 120 MHz by means of hadronically produced secondary elec- 
trons (assuming B = 1 [iG). As expected, weak shocks with Mach 
numbers I < M < 2 almost disappear in this representation due 
to their small acceleration efficiency. This eflfect is amplified by 
considering only energetic CR protons with q > 30. Comparing 
the Mach numbers weighted by ^cr to those weighted by ^diss un- 
covers the three-dimensional picture of these shock surfaces. The 
powerful (in terms of energy dissipation rate) but weak (in terms 
of shock strength) internal shock waves are surrounded by shock 
surfaces of successively increasing Mach numbers that are only re- 
vealed in these projections if one disfavours these internal shocks 
in the weighting function. 

What are the consequences for the CR proton pressure within 
galaxy clusters? Similar to the thermal pressure, it peaks in the cen- 
tre and falls of with radius. The CR pressure is additionally en- 
hanced at strong shock waves, that have a modulating eflfect on the 
underlying CR pressure distribution. The latter can be seen by look- 
ing at the strong tangential shock wave at r ^ 600 kpc south- 
wards from the cluster centre in the mass weighted CR proton pres- 
sure map (top left panel in Fig. 6). The CR pressure peaks roughly 
at PcR - 10~^^erg cm~^ /i^q at the cluster centre. Even more reveal- 
ing is the mass-weighted CR proton pressure relative to the total 
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Figure 5. The environment of a large post-merging galaxy cluster (M ^ 10^^/z~^ M©) in our simulation with radiative gas physics and star formation including 
CRs from structure formation shocks only (model S2). Shown are the line-of-sight averaged density (top left side), the mass weighted temperature (top right 
side), the Mach number of shocks weighted by the energy dissipation rate in colour (while the brightness displays the logarithm of the dissipation rate, bottom 
left side), and the Mach number of shocks weighted by the energy injection rate of CR protons in colour (while the brightness displays the logarithm of the 
CR proton energy injection rate, bottom right side). For better comparison, we used the same colour and brightness scale in the bottom plots. Only CR protons 
with a dimensionless momentum q = fij > 30 have been considered for calculating the CR energy density since lower energetic CR protons are not detectable 
at radio frequencies y > 120 MHz by means of hadronically produced secondary electrons. 



pressure Xcr = Pcr/Piou where Aot = ^cr + Ah (top right panel 
in Fig. 6). The relative CR pressure Xcr acquires comparatively 
high values within the WHIM that are hydrodynamically impor- 
tant, their importance decreases (on average) as we move inwards 
due to a combination of the following reasons: (1) weak central 
flow shocks are inefficient in accelerating CRs (e.g., Paper I) and 
(2) adiabatic compression of a composite of CRs and thermal gas 
disfavours the CR pressure relative to the thermal pressure due to 
the softer equation of state of CRs. Within each individual galaxy as 
well as within the cluster centre, the CR pressure reaches equiparti- 



tion or dominates the thermal pressure as can be seen by the numer- 
ous yellow points sprinkled over the map, each corresponding to a 
galaxy. This is due to the long CR cooling time scales compared to 
those of the thermal gas, an effect that diminishes the thermal gas 
pressure relative to that of CRs (Paper I). 

It is very instructive to compare the CR proton to the CR elec- 
tron pressure since protons and electrons are subject to different 
cooling mechanisms due to their large mass diff'erence. The CR 
proton cooling timescale is generally larger than that of CR elec- 
trons such that protons accumulate within the ICM on a Hubble 
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Figure 6. CR proton and electron pressure maps in the super-cluster environment of a large post-merging galaxy cluster (M ^ 10 /z~ M©) in our radiative 
simulation (model S2). The mass weighted CR pressure (top left side) is contrasted with the mass weighted CR pressure relative to the total gas pressure 
(top right side). Since the CR proton pressure decreases less steeply than the gas pressure as a function of cluster radius this results in an increasing relative 
CR pressure profile towards the periphery. In the bottom panels, we show the mass weighted CR electron pressure relative to the total pressure for primary, 
shock-accelerated electrons (left side) and secondary electrons resulting from hadronic CR proton interactions (right side). The CR electron pressure derives 
from the respective equilibrium distribution functions that balance acceleration and cooling processes. 



timescale and maintain a comparatively smooth distribution over 
the cluster volume (top panels of Fig. 6). This implies that the CR 
proton pressure traces the time integrated non-equilibrium activi- 
ties of a cluster and is only modulated by recent dynamical activ- 
ities (see also Paper I, for average values of the relative CR en- 
ergy in different dynamical cluster environments). In contrast, the 
pressure of primary CR electrons resembles the current dynami- 
cal, non-equilibrium activity of the forming structure and results 
in an inhomogeneous and aspherical spatial distribution. To under- 
pin this argument, in the bottom panels of Fig. 6, we show the 
mass-weighted CR electron pressure relative to the total pressure 



^CRe = ^CRe/Aot- On the left side, we show the relative pressure 
of primary, shock- accelerated electrons while the relative pressure 
of secondary electrons resulting from hadronic CR proton interac- 
tions is shown on the right side. The CR electron pressure derives 
from the respective equilibrium distribution functions that balance 
acceleration and cooling processes as laid out in Appendices A3. 2 
and A4.2. Note that the colour scale for the relative pressure of pri- 
mary CR electrons spans exactly two orders of magnitude (like in 
the case of CR protons), peaking in the dilute WHIM at roughly 
3 per cent rather than 30 per cent as in the case of CR protons. 
The relative CR electron pressure XcRe,prim decreases towards clus- 
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ters and groups due to larger Coulomb losses and smaller shock 
acceleration efficiencies within collapsed objects as in the case of 
CR protons. Interestingly, XcRe,prim within galaxies is suppressed by 
roughly one order of magnitude with respect to the ambient inter- 
galactic medium in which the galaxy resides due to large Coulomb 
losses. This is quite diff'erent from CR protons that acquire equipar- 
tition with the thermal gas inside galaxies. In general, the spatial 
variations of XcRe,prim are larger than in the case of protons, showing 
that the CR electron pressure indeed reflects the active dynamical 
structure formation activities mediated by shock waves. In contrast, 
the mass weighted relative pressure of secondary CR electrons is 
shown on a colour scale that spans seven orders of magnitude, due 
to the high dynamic range of this quantity. The CR electron pres- 
sure is proportional to the number densities of CR protons and of 
the gas PcRcsec ^ ^cr'^n^ causing XcRe,sec to peak towards dens- 
est structures and thus filling in the diminishing primary CR elec- 
tron pressure inside dense structures such as galaxies. As a word of 
caution, we do not account for the re- acceleration of CR electrons 
e.g. via resonant pitch angle scattering by compressible magneto- 
hydrodynamical (MHD) modes neither do we account for a pre- 
viously injected and aged electron population which could change 
the presented picture. Further work is required to elucidate these 
electron components in simulations. 

3.2 Radio synchrotron emission 

3. 2. 1 Projected radio maps 

Figure 7 shows the large-scale "radio web"of the super-cluster re- 
gion of our Coma-like cluster that experienced a recent merger. The 
radio emission is computed assuming a simple scaling model for 
the magnetic field of 



where the central magnetic energy density ssfi and as are free pa- 
rameters in our model, and £:th,o denotes the thermal energy density 
at the cluster centre. It is motivated by turbulent dynamo theories 
of the growth of magnetic field strength that are typically saturat- 
ing at a level proportional to the turbulent energy density or the 
thermal energy density (e.g., Subramanian 2003; Schekochihin & 
Cowley 2006). This allows us to explore the unknown behaviour of 
the large scale magnetic field parametrically (cf . Appendix B 1 for 
more discussion). Our standard model (also adopted in Fig. 7) as- 
sumes a central magnetic field strength of 5o = 10 \xG and of^ = 0.5 
which implies a constant ratio of thermal-to-magnetic pressure of 
20 in our Coma-like cluster. 

In the top panels of Fig. 7, we separately show the synchrotron 
emission of primary CR electrons that were accelerated directly at 
structure formation shocks as well as the radio emission of sec- 
ondary CR electrons that results from hadronic CR proton interac- 
tions with ambient gas protons. The combined radio synchrotron 
emission (shown in the bottom left panel) shows that the morpho- 
logically smooth secondary component dominates the radio emis- 
sion of the central cluster regions. In contrast, the irregularly shaped 
primary radio relic emission dominates in the cluster periphery and 
the super-cluster region that is believed to host the warm-hot inter- 
galactic medium (WHIM). These observations are supported by the 
radio spectral index map between 150 MHz and 1.4 GHz (bottom 
right panel) which shows larger variations in the peripheral cluster 
regions. These are caused by projecting the radio emission from 
inhomogeneously accelerated primary CR electrons and reflect the 



strong variation of the Mach numbers of structure formation shocks 
at the outer cluster regions. Based on these findings, we put for- 
ward a new unified scheme for the generation of giant radio halos 
as well as radio mini-halos. The giant radio halo emission in merg- 
ing clusters shows a transition from the secondary radio emission 
in the centre to the dominant primary emission component at the 
outer parts of radio halos. Gravitational energy, that is associated 
with the merger, is virialised by a morphologically complex net- 
work of strong shock waves in the cluster outskirts. This induces an 
irregular radio 'gischt' emission in the cluster periphery that repre- 
sents radio synchrotron radiation emitted from shock-accelerated 
electrons. Our simulated radio emission maps of relaxed cool core 
clusters show a significantly reduced level of this primary emission 
component such that the diffuse radio emission in these systems is 
solely determined by the secondary radio emission, producing a ra- 
dio mini-halo. Note that our simple magnetic model does not take 
into account the adiabatic compression of magnetic fields during 
the formation of a cool core. This effect should furthermore de- 
crease the emission size of radio mini-halos making it comparable 
to the cool core region. 

Closer inspection of the primary radio emission map (top left 
panel in Fig. 7) shows a bright radio relic on the lower right with 
respect to the cluster centre at a distance of r ^ 0.6/z~^ Mpc. This 
is caused by an outgoing merger shock wave that steepens as it 
reaches the shallower peripheral cluster potential. Further outwards 
at a distance of 2 and 3 h~^ Mpc to the cluster centre, there is a class 
of tangentially curved radio relics visible in orange and red. These 
are uniquely associated with strong shock waves as can be inferred 
by comparing the primary radio emission to the dissipated energy 
at shock waves (shown as brightness in the bottom left panel of 
Fig. 5). The statistical study of the radio emission at different fre- 
quencies and Faraday rotation of these objects will enable us to 
investigate non-equilibrium processes of virialisation including the 
acceleration of cosmic rays, the growth of magnetic fields, and ki- 
netic energy in bulk motions that are expected to source turbulence 
in clusters. Comparing the emission level of the projected radio 
surface brightness maps to the LOFAR point source sensitivity of 
0.25 mJy / (arcmin hour) at y = 120 MHz shows that moderately 
long exposures of super cluster regions have the potential to de- 
tect the large-scale "radio web" and to study the magnetic field on 
Mpc- scales that is woven into the web. Ongoing work that includes 
simulated mock observations for radio interferometers studies as- 
sociated questions in greater detail (Battaglia et al. in prep.). 

Figure 8 shows the dependence of the primary and secondary 
radio synchrotron emission on the observing frequency, the model 
for the magnetic field, as well as the type of simulated gas physics 
(radiative versus non-radiative). The top panels show the high- and 
low-frequency radio emission (1.4 GHz and 15 MHz) in our radia- 
tive simulation (model S2) using our standard model for the mag- 
netic field. This demonstrates the potential of low-frequency ra- 
dio arrays in studying non-thermal properties of the inter-galactic 
medium especially since the associated radio spectrum is steeper 
compared to that of the Galactic foreground emission. This will al- 
low us to address questions such as the existence and properties 
of the WHIM and the existence and origin of large scale mag- 
netic fields. The bottom left panel shows the total radio emission 
at 150 MHz in our non-radiative simulation (model SI) assuming 
our standard parameters for the magnetic field, and should be com- 
pared to the same panel in Fig. 7. The level of the primary radio 
emission in the cluster periphery and the super-cluster region is re- 
duced in the non-radiative simulation (model SI) compared to the 
radiative case (model S2). Some relics in the bottom panel (model 
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Figure 7. The large-scale "radio web" at 150 MHz of the super-cluster region of our Coma-like cluster that experienced a recent merger (g72a) in our radiative 
simulation (model S2). We show the synchrotron emission of primary CR electrons that were accelerated directly at structure formation shocks (top left side) as 
well as the radio emission of secondary CR electrons that results from hadronic CR proton interactions with ambient gas protons (top right side). The bottom 
left panel shows the giant radio halo emission of this cluster that is characterised in the centre by the regular smooth morphology of the secondary radio 
emission. At larger radii, we observe a transition to the the irregularly shaped primary radio "gischt" emission with a prominent radio relic to the lower right 
of the cluster. The radio spectral index between 150 MHz and 1.4 GHz (bottom right panel) shows larger variations in the peripheral cluster regions. These 
are caused by projecting the radio emission from inhomogeneously accelerated primary CR electrons and reflect the strong variation of the Mach numbers of 
structure formation shocks at the outer cluster regions. 



SI) that are at distances of > 2h~^ Mpc to the cluster centre even 
show radio emission at the level that is comparable to that in the 
panel above (model S2) despite the lower frequency that should 
provide a flux level that is increased by an order of magnitude as- 
suming Qfy ^ 1. The reason for this stems from the larger shock 
strength (higher Mach numbers) of characteristic shocks that dis- 
sipate gravitational energy into thermal energy in radiative simu- 
lations (Paper I). The enhanced acceleration efficiency of CRs at 
stronger shocks leads to the increased primary radio emission in 
radiative simulations compared to the non-radiative case. The bot- 



tom right panel again shows the total radio emission at 150 MHz in 
our non-radiative simulation, however, with a shallower magnetic 
decline, aB - 0.25 which results in sb ^ s^^^ . Although the radial 
decline of this model for the magnetic field might be almost too 
shallow, it serves for illustrative purposes demonstrating that low- 
frequency radio arrays in combination with high-resolution simula- 
tions can tightly constrain the large scale behaviour of the magnetic 
field. 
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Figure 8. Dependence of the primary and secondary radio synchrotron emission on the observing frequency, the model for the magnetic field, as well as 
the type of simulated gas physics (radiative versus non-radiative). The top panels show the synchrotron radiation at 1.4 GHz and 15 MHz in our radiative 
simulation (S2) assuming a simple scaling of the magnetic energy density with the thermal energy density, sb oc st\i. This demonstrates the potential of low- 
frequency radio arrays in studying non-thermal properties of the inter-galactic medium. The bottom left panel shows the total radio emission at 150 MHz in 
our non-radiative simulation (SI) with the same magnetic model. The bottom right panel again shows the total radio emission at 150 MHz in our non-radiative 
simulation, however, with a shallower magnetic decline, sb oc s^.^. For convenience, the colour scale is the same in all panels such that the emission at 15 
MHz in the central cluster region is highly saturated. 



3.2.2 Radio emission profiles 

Radio synchrotron profiles allow us to confirm and quantify the 
proposed unified scheme for the cluster radio halo emission. 

Primary versus secondary radio emission: The left side of 
Fig. 9 compares the synchrotron emission of primary CR electrons 
that were accelerated directly at structure formation shocks with 
that of secondary CR electrons that result from hadronic CR pro- 
ton interactions with ambient gas protons. In cool core clusters, 
the azimuthally averaged secondary radio emission dominates the 
primary emission component for radii r < 3/?vir- The smooth sec- 



ondary component typically falls of at a radius r ^ 0.2/?vir which 
resembles the characteristics of observed radio mini halos as ob- 
served e.g. in the Perseus cluster (Pedlar et al. 1990) or RX J1347.5- 
1 145 (Gitti et al. 2007). Our model predicts diflfuse secondary radio 
emission in virtually every cool core cluster. 

Interestingly, our post-merging cluster g72a shows a transi- 
tion from the secondary to the primary radio emission component 
towards the outer cluster regions triggered by the dynamical merger 
activity with strong shock waves traversing the cluster in order to 
thermalize the gas. These shock waves steepen as they break on 
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Primary vs. secondary radio emission: Influence ofsim. physics and magnetic models: 




Figure 9. Azimuthally averaged radio synchrotron brightness profiles at y = 1.4 GHz. The left side compares the synchrotron emission of primary CR electrons 
that were accelerated directly at structure formation shocks (dotted lines) with that of secondary CR electrons that result from hadronic CR proton interactions 
with ambient gas protons (dashed line) while the soHd line shows the sum of both emission components. The upper panel compares a post-merging cluster 
(g72a) with a cool core (CC) cluster (g51) of the same mass and the bottom panel compare the radio emission of two differently sized CC clusters (g8a versus 
g676). The right side shows the influence of different simulation and magnetic field models on the total synchrotron emission profiles for the post-merging 
cluster g72a (upper panel) and the CC cluster g8a (bottom panel). The standard model uses radiative gas physics (S2), has a central magnetic field strength of 
5o = 10 ia.G, and scales as oc s^^^^ with ag = 0.5. The other lines result from varying each of these assumptions separately leaving the others unchanged. 
The radio emission in our non-radiative simulation (SI, dotted lines) declines faster. The emission profiles for a lower central magnetic field strength (dashed 
lines) and with a weaker magnetic decline (dash-dotted lines) illustrate the uncertainty in the simulated radio profiles. 



the shallower peripheral cluster potential. The associated increase 
in the primary radio emission comes hand in hand with a stronger 
variation of the radio spectral index towards the cluster periph- 
ery (cf. Fig. 7). In the particular case of g72a, the resulting ra- 
dio halo profile reaches out to r ^ 0.45/?vir which corresponds to 
the observed size of /?max - 1 Mpc/Zzyo of the Coma radio halo at 
y = 1.4 GHz (Deiss et al. 1997; Reiprich & Bohringer 2002). We 
verified that the transition from the secondary to the primary radio 
emission in our simulated giant radio halos is independent on the 
chosen projection and a generic prediction for merging clusters. 
Our simple scaling model for the magnetic field of Eqn. (5) does 
clearly not include non-equilibrium effects related to the growth of 
the magnetic field. The enhancement of the magnetic field strength 
through turbulent dynamo processes will saturate on a level which 
is determined by the strength of the magnetic back-reaction (e.g., 
Subramanian 2003) and is typically a fraction of the turbulent en- 
ergy density. Thus, in a real cluster, the strong shocks at the cluster 
periphery are expected to drive turbulence and strong shear motions 
which should in turn lead to a stronger magnetic field amplification. 
Our adopted scaling of the magnetic field with the thermal energy 



density might partially neglect these effects and should somewhat 
underestimate the peripheral radio synchrotron emission. 

Influence of magnetic parametrisation on radio emission: 

The right side of Fig. 9 shows the influence of different simula- 
tions and magnetic field models on the total synchrotron emission 
profiles for our massive post-merging cluster and cool core clus- 
ter. The emission profiles for lower central magnetic field strength 
(dashed lines, Bq = 3 \iG) shows a small decrease of the central ra- 
dio emission by a factor of two while it is considerably suppressed 
by an order of magnitude towards larger radii. This is due to the 
two distinctive regimes of a synchrotron emitting equilibrium dis- 
tribution of CR electrons (cf. Fig. 3). For field strengths B > 3 \iG, 
the synchrotron flux is almost insensitive to the field strength while 
it scales as jy oc for weaker fields B <3 ]iG which are present 
at larger cluster radii. The emission profile for a weaker magnetic 
decline (dash-dotted lines, as = 0.25) is more extended than our 
standard model, as expected. Current cosmological MHD SPH sim- 
ulations (e.g., Dolag et al. 2001) may not sufficiently resolve small 
scale turbulent dynamo processes and large shear motions that are 
thought to amplify the magnetic field in the coarsely sampled super- 
cluster regions beyond the accretion shocks. The model with the 
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Figure 10. Average profiles for the spherically averaged CR fraction /cr = «CR/«th of our sample of all 14 clusters at redshift z = 0. The error bars represent 
the standard deviation from the sample mean. The increase of /cr with radius is independent of the modelled CR physics (SI - solid, S2 - dashed, S3 - 
dash-dotted) and continues beyond the virial radius suggesting that this is a generic property of CR physics. The increase of /cr in our radiative models for 
small radii r < 0. l/?vir is due to the short cooling timescale of the gas compared to that of CR protons (cf. Paper I). Note the different axes scales in both panels. 



weaker magnetic decline is an attempt to parametrise those un- 
certainties. Note that despite the uncertainties in the parametrisa- 
tion of the magnetic field and thus the overall radio emission, the 
conclusions with respect to the different emission components (pri- 
mary versus secondary) and their emission characteristics remain 
unchanged. 

Influence of simulated physics on radio emission: The radio 
emission in our non-radiative simulations (SI, dotted lines) is much 
smoother and declines faster since there is only a very weak transi- 
tion to the primary component due to the weaker shocks in the clus- 
ter periphery compared to the radiative simulations (S2). It is sur- 
prising that the central radio emission in simulation models SI and 
S2 almost coincide despite the large difference of the central CR 
fraction /cr = ^cR/^th in both models (cf. Fig. 10). At r = 0.02/?vir 
the CR fraction is two orders of magnitude larger in our radiative 
simulations compared to our non-radiative simulations. This can be 
understood by the self-regulated nature of CR feedback. The sec- 
ondary synchrotron emission scales as jy oc /icR^gas ^ /cR^gas 
glecting the weak additional density dependence through the mag- 
netic field in the synchrotron regime (cf. Fig. 3). The lower gas 
density in the radiative simulations (cf. Figs. 3 and 5 in Paper I) 
almost exactly balances this difference of the CR fraction such that 
the resulting secondary synchrotron emission level in the centre re- 
mains only slightly modified. This is due to a combination of the 
following reasons. (1) The CR cooling timescales due to Coulomb 
and hadronic interactions of CRs, Tpp,coui ^ ^gas' almost an or- 
der of magnitude larger in our non-radiative simulations compared 
to our radiative case owing to the central density difference. (2) 
A second sub-dominant effect is the reduced depletion of the CR 
pressure in our radiative simulations due to adiabatic compression 
of our composite of CRs and thermal gas which disfavours the CR 
pressure relative to the thermal pressure. 

Part of this density difference is reinforced in contemporary 
cosmological radiative simulations that do not include feedback 
from AGN. This leads to the well-known over-cooling problem 
which results in an overproduction of the amount of stars, enhanced 
central gas densities, and too small central temperatures compared 
to X-ray observations. The density enhancement at the very centre 



and the associated star formation take place at the expense of the 
surrounding ICM which ends up being less dense compared to its 
initial stage before cooling set in. This hypothetical initial stage is 
realised by our non-radiative simulations that does not take into ac- 
count radiative cooling. We show that the secondary CR emission 
(radio synchrotron, inverse Compton, and pion decay induced y-ray 
emission) within the framework of our CR model is almost inde- 
pendent of those short-comings in the central cluster regions. The 
difference of the radio emission at larger radii between our models 
SI and S2 however is a robust finding and primarily caused by the 
difference of the primary radio emission. This difference is due to 
the on average stronger shock waves that lead to more efficient CR 
electron acceleration in our radiative simulations. 



3.2.3 Discussion of synchrotron polarisation 

Synchrotron emission from primary accelerated electrons is po- 
larised due to a combination of two effects. (1) Shock compression 
aligns unordered magnetic fields with the shock plane (Ensslin et al. 
1998). Furthermore, shearing motions induced by oblique shocks 
stretch these field lines which leads to a larger magnetic coher- 
ence length of two-dimensional field configurations (Schekochihin 
& Cowley 2006). If the synchrotron emitting structure that is ener- 
gised by the formation shock wave is seen at some angle between 
the line-of- sight and the normal of the shock front, the magnetic 
field structure projected onto the plane of sky shows a preferential 
direction which implies a preferred intrinsic synchrotron polarisa- 
tion. (2) The combination of the localised acceleration site of CR 
electrons at these shock fronts and the short synchrotron cooling 
times (cf. Fig. Al) leads to a small synchrotron emission volume. 
Thus, these peripheral radio relics are expected to show a preferred 
synchrotron polarisation with the magnetic field aligned with the 
shock surface (as observed e.g. in Abell 3667 by Rottgering et al. 
1997). Superposing many causally unconnected radio relics in pro- 
jection leads to a decrease of the degree of polarisation. 

Hadronically induced synchrotron emission of the smooth ra- 
dio halo is virtually unpolarised assuming statistically isotropic dis- 
tribution of magnetic field vectors without a preferred direction. 
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Thermal X-ray emission: Pion decay y-ray emission (Ey > 100 MeV): 
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Figure 11. The top panels compare the thermal X-ray emission and the y-ray emission resulting from hadronic CR interaction with ambient gas protons of 
the super-cluster region of our Coma-like cluster in our radiative simulation (model S2, cluster g72a). The hadronic y-ray emission shows a shallow decline 
with radius due to the rising CR-to-thermal number density profile. The bottom panels show the inverse Compton emission from primary and secondary CR 
electrons in the hard X-ray (left side) as well as the y-ray band (right side). The primary CR electrons dominate the emission signal on large scales. Comparing 
the y-ray emission components (right panels) shows that the pion decay y-rays exceed the total IC emission at energies > 100 MeV. 



The large emission volume is filled with magnetised plasma that 
causes the plane of polarisation to Faraday rotate. Hadronically 
generated CR electrons fill the same cluster volume. Thus, each 
radio emitting volume element along the line-of- sight, that is sep- 
arated by more than the magnetic correlation length or the Fara- 
day depth, if shorter, radiates causally unconnected intrinsically 
polarised emission that averages out to a net unpolarised emission, 
e.g. radio halos are Faraday depolarised. 

Combining these considerations with the previously devel- 
oped model for the radio halo emission implies a transition from 
the virtually unpolarised radio halo emission at small impact pa- 
rameters to a small degree of polarisation at the halo periphery. 



characterised by the dominating primary emission there. In order to 
detect this polarisation one might be forced to go out to large im- 
pact parameters with a small resulting synchrotron surface bright- 
ness where the emission is dominated by very few relics along the 
line-of- sight. Owing to the diff'erent injection timescales of primary 
and secondary CR electrons, we conclude that the secondary halo 
emission traces the time integrated non-equilibrium activities of a 
cluster and is modulated by the recent dynamical activities (Paper 
I). In contrast, the polarised radio relic emission resembles the cur- 
rent dynamical, non-equilibrium activity of a forming structure and 
results in an inhomogeneous and aspherical spatial distribution. 
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Different y-ray emission processes: 



Comparison of different clusters: 




Figure 12. Azimuthally averaged profiles of the inverse Compton (IC) surface brightness for energies > 10 keV (upper panels) and y-ray emission profiles 
for Ey > 100 MeV (bottom panels). Using the simulation of the post-merging cluster g72a, the left side compares the IC emission of primary CR electrons 
that were accelerated directly at structure formation shocks (dotted lines) with that of secondary CR electrons that result from hadronic CR proton interactions 
with ambient gas protons (dashed lines) while the solid lines shows the sum of both emission components. The dominant emission component at energies 
Ey > 100 MeV is the pion decay induced y-ray emission (dash-dotted lines). The right side compares the total hard X-ray and y-ray emission, respectively, 
for for the different clusters g8a (dotted lines), g72a (solid lines), g51 (dashed lines), and g676 (dashed-dotted lines). While the hard X-ray/y-ray emission 
clearly scales with the cluster mass, the dynamical state of the cluster is equally important and can even reverse the mass trend as can be inferred from the IC 
emission for Ey > \0 keV comparing the top two clusters in the figure legend. 



3.3 Hard X-ray and y-ray emission 

3. 3. 1 Projected X-ray and y-ray maps 

Different non-thermal y-ray emission processes are compared to 
the thermal X-ray emission in Fig. 11. The top panels compare 
the morphology of the thermal X-ray surface brightness (for de- 
tails of the projection, cf. Paper I) to that of the y-ray emission re- 
sulting from hadronic CR interactions with ambient gas protons of 
our Coma cluster region in our radiative simulation. Although they 
resemble each other, the pion decay induced y-ray emission de- 
clines slower with radius and makes the y-ray halo more extended. 
The bremsstrahlung emission of the thermal gas scales as (ne- 
glecting the radial dependence of the cooling function) while the 
hadronically induced y-ray emission scales as Wgas^cR- The discrep- 
ancy between the radial behaviour of the thermal and non-thermal 
emission can be understood by looking at the right side of Fig. 10 
which shows average profiles of the spherically averaged CR frac- 
tion /cR = ficR/nth of our sample of all clusters. There is an in- 
crease of fcR with radius that is independent of the modelled CR 
physics and continues beyond the virial radius. This generic pre- 
diction of CR physics is due to the more efficient CR acceleration 



at the peripheral strong accretion shocks compared to weak central 
flow shocks. 

The bottom panels of Fig. 1 1 show the inverse Compton emis- 
sion from primary and secondary CR electrons in the hard X-ray 
(left side) as well as the y-ray band (right side). The primary CR 
electrons dominate the hard X-ray emission signal on large scales. 
The primary IC emission directly reflects the inhomogeneous viri- 
alisation process that manifests itself through a filigree web spun by 
shocks (cf. Fig. 5). In principle, IC emission is the cleanest way of 
probing structure formation shock waves since the inverse Comp- 
ton emission is not weighted by the magnetic energy density as it 
is the case for synchrotron emission. Visually comparing the right 
panels of Fig. 1 1 implies that the pion decay induced y-ray emis- 
sion dominates the total IC emission in the energy range that is of 
interest to GLAST. 



3. 3. 2 Hard X- ray IC and y- ray emission profiles 

Figure 12 shows azimuthally averaged profiles of the IC surface 
brightness for energies Ey > 10 keV (upper panels) and y-ray emis- 
sion profiles for Ey > 100 MeV (bottom panels). Using the simu- 
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lation of the post-merging cluster g72a, the left side compares the 
IC emission of primary CR electrons that were accelerated directly 
at structure formation shocks with that of secondary CR electrons 
that result from hadronic CR proton interactions with ambient gas 
protons. The primary and especially the secondary IC emission is 
generally more inhomogeneous compared to both radio emission 
components. Secondly, the IC emission declines less steeply with 
radius and shows almost a power-law profile compared with the p- 
profile in the case of synchrotron emission. Both efi'ects are due to 
the weighting by the magnetic field which in the latter case tends to 
smooth out the inhomogeneous CR electron distribution and causes 
a steeper decrease of the radio synchrotron emission (at least in our 
model of the magnetic field). 

At energies Ey > 100 MeV, the pion decay induced y-ray 
emission is the dominant emission component everywhere except 
in the peripheral regions of the post- merging cluster g72a where the 
primary IC emission is of similar strength (cf. bottom left panel of 
Fig. 12). As a word of caution it should be added, that our simula- 
tion assumes a CR proton spectral index of ofp = 2.3. Lowering this 
value would result in a larger secondary IC component (cf. Fig. 4). 
Future work is required to study such a scenario which might find 
applications in the outer parts of clusters. However, the pion de- 
cay component at Ey > 100 MeV should be robust with respect 
to variations of since it samples the pion bump that is sensi- 
tive to threshold eff'ects of the hadronic cross-section but not to the 
spectrum of the parent CR distribution. Despite the fact that the 
secondary IC emission still dominates the primary IC emission at 
Ey > 100 MeV, the primary IC component is increased by a fac- 
tor of two compared to the hard X-ray emission at Ey > 10 keV. 
We conclude that the mean IC spectral index is obviously smaller 
than that of the secondary emission, ofy = 1.15. This can be un- 
derstood by combining the facts that a superposition of diff'erent 
power-law spectra produces a concave spectrum and that IC emis- 
sion at Ey > 100 MeV results from CR electrons with a Lorentz 
factor 7 > 3 X 10^ according to Eqn. (2). This again stresses the 
importance of correctly modelling the peripheral regions of a clus- 
ter since they show predominantly the conditions for strong shock 
waves that are able to accelerate such flat CR electron populations. 

We move on to compare the hard X-ray and y-ray emission for 
clusters of diff'erent masses and dynamical states. While the hard X- 
ray/y-ray emission clearly scales with the cluster mass, the dynam- 
ical state of the cluster is equally important and can even reverse 
the mass trend (cf. the right side of Fig. 12). The reason for this 
lies in the enhanced CR pressure in merging clusters (Paper I) as 
well as in the primary IC emission that sensitively traces current 
non-equilibrium or merging activities of clusters. 

Figures 14 and 13 study the influence of diff'erent simulated 
physics on y-ray emission at Ey > 100 MeV. The y-ray emission 
is more inhomogeneous in radiative simulations compared to non- 
radiative simulations. This is due to the CR pressure equipartition 
within each galaxy. The overall y-ray luminosity, however, is very 
similar as can be inferred from the azimuthally averaged emission 
profiles. This confirms our finding for the synchrotron emission 
in Fig. 9 and confirms our explanation that this is indeed a self- 
regulated efl'ect of CR feedback and not biased due to the magnetic 
weighting of the synchrotron emission. The complete CR model 
that also accounts for CR from SNe (dashed lines, model S3) shows 
a slightly enhanced y-ray emission level since we opened up a new 
acceleration channel compared to our standard model S2. 



Influence of simulated physics on y-ray emission: 
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Figure 13. Influence of simulated physics on azimuthally averaged y-ray 
emission profiles for Ey > 100 MeV. The top panel shows the emission 
profile of the post-merging cluster g72a while the bottom panel that of the 
large cool core cluster g8a. Shown are our non-radiative simulations (dot- 
ted lines, model SI), our radiative simulations with CR protons injected at 
structure formation shocks (solid lines, model S2), and the complete CR 
model where we additionally take CR acceleration at SNe shocks into ac- 
count (dashed lines, model S3). Despite the diff'erent physics in these simu- 
lations, the y-ray emission level is similar. 



3.4 Correlations between thermal and non-thermal emission 

Compressing intrinsically non-spherical emission purely into sur- 
face brightness profiles causes loss of information and might yield 
biased results. In addition to emission maps and surface brightness 
profiles, we complete our analysis using pixel-to-pixel correlations 
of the thermal X-ray surface brightness with non- thermal cluster 
emission processes. To this end we compare these correlations of 
our post-merging cluster simulation g72 (Fig. 15) to the cool core 
cluster simulation g8 (Fig. 16). Each of these figures shows the 
correlation space density of the radio surface brightness (top pan- 
els) and the y-ray surface brightness for Ey > 100 MeV (bottom 
panels), as well as that of the hadronically induced non- thermal 
emission (left side, red colour scale) and the non-thermal emission 
of primary CR electrons that were accelerated directly at structure 
formation shocks (right side, blue colour scale). 

While the hadronically induced non-thermal emission is 
tightly correlated with the thermal bremsstrahlung emission, the 
correlation is much weaker and the scatter is increased in the case 
of primary non-thermal emission where structures in the correlation 
space density correspond to individual structure formation shock 
waves. We can see preferably tangential shocks that are charac- 
terised by a varying non-thermal emission for a constant X-ray sur- 
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Radiative simulation, y-rays (Ey > 100 MeV): 




Non-radiative simulation, y-rays (Ey > 100 MeV): 




Figure 14. Influence of simulated physics on pion decay y-ray emission maps for Ey > 100 MeV of the large cool core cluster g8a. The peaked emission 
of each galaxy in the radiative simulation (left panel, model S2) contrasts the smoother and fluffier y-ray emission of our non-radiative simulation with CR 
protons injected at structure formation shocks (right side, model SI). 



face brightness and to a smaller extend radial shocks where the role 
of thermal and non-thermal emission is interchanged. 



3.4.1 Correlations of the synchrotron emission 

Closer inspection of the secondary synchrotron emission (top left 
panel in Fig. 15) shows flattening of the correlation 



(6) 



where the power-law index changes from a - 1.7 to a = 1.3 



above Sxfi = 3x10 ^ erg cm 



hlr., with a normalisation Sy 



0.057 mJy arcmin" /z^^. This is due to a combination of the fol- 
lowing two eff'ects. (1) Towards higher X-ray surface brightness, 
the radio emission experiences a transition from IC to synchrotron 
regime with a weaker dependence on magnetic field that goes along 
with weaker density dependence. (2) The merger displaces the cen- 
tral cool core and disturbs the ICM by means of merging shock 
waves that dissipate the gravitational binding energy associated 
with the merger. This yields an increased CR proton pressure and 
number density relative to the thermal gas within the central X- 
ray luminous regions and causes the change in the correlation. 
Since the thermal energy distribution is equally eff'ected by these 
re-distribution of energy, this merger induced eff'ect can be further- 
more amplified by the transition from the IC to synchrotron regime 
as explained above. The change of the power-law index of the cor- 
relation is more pronounced for our post-merging cluster simula- 
tion g72 compared to our cool core cluster g8 suggesting the im- 
portance of the second eff'ect. Note that the explicit values of the 
correlation, in particular 5" yo» depend on the assumed model for the 
magnetic field. There is a second branch in the correlation of the 
hadronically induced non-thermal emission visible that is due to 
the other smaller cluster forming in that simulation. 

The correlation of the primary synchrotron emission (top right 
side) shows a large scatter especially at large surface brightness 



such that the correlation of the total synchrotron brightness is ex- 
pected to broaden and to become flatter towards lower surface 
brightness to yield a roughly linear correlation. The exact realisa- 
tion of the correlation at a surface brightness that is substantially 
supported by primary synchrotron emission sensitively depends on 
the mass ratio, geometry, and the advanced state of the merger. 
Thus our simulations can only provide quantitative predictions for 
the statistical behaviour rather than deterministic predictions for 
the correlations. Our correlations are strikingly similar to the ones 
found in observed radio halos (cf. Govoni et al. 2001). They find 
a linear relation between the radio and X-ray surface brightness 
that is tight at high surface brightness while it broadens and flattens 
towards dimmer brightness. This behaviour is one of the strongest 
arguments in favour of our new model for radio halos that should be 
dominated in the centre by secondary emission with a transition to 
the primary synchrotron emission in the cluster periphery. In accor- 
dance with our findings, their radio halo emission is slightly more 
extended compared to the thermal X-ray emission. Varying spec- 
tral index distributions preferably in the cluster periphery (Feretti 
et al. 2004) support this picture. In particular, our model supports 
a strong link between radio halos and cluster mergers for which 
there is a strong evidence in the literature (Feretti et al. 2004, and 
references therein). 



3.4.2 Correlations of the y-ray emission 

The bottom left panels in Figs. 15 and 16 include pion decay y- 
ray emission as well as IC emission from secondary CR electrons. 
These tight correlations are characterised by a sub-linear power- 
law relation. This is due to the shallower decay of the CR number 
density ncK compared to that of the thermal gas leading to an in- 
crease of the CR fraction ncR/n± (cf. Fig. 10). Our post-merging 
cluster g72a shows a small variation of the power-law index of the 
y-ray to X-ray correlation at high surface brightness albeit not as 
pronounced as in the case of radio synchrotron emission. This con- 
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Figure 15. Pixel-to-pixel correlation of the thermal X-ray surface brightness with both the radio surface brightness (top panels) and the y-ray surface brightness 
for Ey > 100 MeV (bottom panels) for our post-merging cluster simulation g72 using our model S2. Shown is the correlation space density of the hadronically 
induced non-thermal emission (left side, red colour scale) and the non-thermal emission of primary CR electrons that were accelerated directly at structure 
formation shocks (right side, blue colour scale). The bottom left panel includes pion decay y-ray emission as well as IC emission from secondary CR electrons. 
While the hadronically induced non-thermal emission is tightly correlated with the thermal bremsstrahlung emission, the correlation is much weaker and the 
scatter is increased in the case of primary non-thermal emission where the structures in the correlation space density correspond to individual structure 
formation shock waves. The line in the top left panel is a fit to the correlation where the slope flattens from 1.7 to 1.3 towards high luminosities. 



firms that both the merger induced boost of the CR pressure and 
the transition from IC to synchrotron regime is responsible for the 
flattening of the correlation between radio and X-ray surface bright- 
ness. The primary IC emission of the cool core cluster simulation 
g8a (lower right panel of Fig. 16) shows a complete absence of any 
correlation. The influence of the merger activity of a cluster on en- 
hancing the non-thermal cluster emission can thus unchallengeably 
be studied with the primary IC y-ray emission (lower right panel 
of Fig. 15). Comparing this primary IC y-ray emission to its coun- 
terpart radio emission indicates that a large part of that correlation 



with the X-ray emission is indeed owed to the merger induced CR 
enhancement and only in parts by the density dependence of the 
magnetic field. 



4 DISCUSSION AND FUTURE WORK 
4.1 Comparison to previous literature 

There have been a series of pioneering papers simulating the non- 
thermal emission from clusters by numerically modelling discre- 
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Figure 16. Same as Fig. 15 but for our cool core cluster simulation g8. Note the absence of any correlation in the bottom right panel showing the primary IC 
surface brightness for Ey > 100 MeV which suggests that the correlation of the primary radio emission (top right panel) is solely due to our assumed model 
of the magnetic field, sb oc s±. 



tised CR energy spectra on top of Eulerian grid-based cosmologi- 
cal simulations (Miniati 2001; Miniati et al. 2001a,b; Miniati 2002, 
2003). In contrast to our approach, these models neglected the hy- 
drodynamic pressure of the CR component, were quite limited in 
their adaptive resolution capability, and they neglected dissipative 
gas physics including radiative cooling, star formation, and super- 
nova feedback. Comparing the non-thermal emission characteris- 
tics of primary CR electrons, hadronically generated secondary CR 
electrons, and pion decay y-rays, we confirm the general picture 
put forward by these authors while we find important diflferences 
on smaller scales especially in cluster cores. Our inhomogeneous, 
peripheral radio relic emission resembles their findings. However, 
the hadronic component of our simulated radio halos is more cen- 
trally concentrated (cf. Miniati et al. 2001b). Our simulations both 



agree that the predicted level of hard X-ray inverse Compton emis- 
sion falls short of the claimed detection in Coma and Perseus albeit 
the discrepancy is more dramatic in our simulations. We confirm 
that the high-energy y-ray emission (Ey > 100 MeV) from clus- 
ter cores is dominated by pion decays while at lower energies, the 
IC emission of secondary CR electrons takes over (Miniati 2003) 
- at least for non-merging clusters. We reproduce their finding that 
the y-ray emission in the virial regions of clusters and beyond in 
super-cluster regions stems from IC emission of primary shock ac- 
celerated electrons. Contrarily to these authors, we find that the sur- 
face brightness of this emission component remains sub-dominant 
in projection compared to the hadronically induced emission com- 
ponents in the cluster core and that the pion decay completely dom- 
inates the high-energy y-ray emission of clusters. We note that our 
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y-ray fluxes from clusters are typically a factor of two smaller than 
the estimates given in Miniati et al. (2001a). As worked out in Pa- 
per III, this has important implications for the number of detectable 
y-ray clusters by GLAST. 

All the discrepancies can be understood by two main eff'ects 
that lead to an overestimation of the CR pressure inside the clus- 
ters simulated by Miniati et al. (2001a) and thus overproduced the 
resulting non-thermal emission: (1) Miniati et al. (2000) identified 
shocks with Mach numbers in the range 4 < At < 5 as the most 
important in thermalizing the plasma. In contrast, Ryu et al. (2003) 
and Pfrommer et al. (2006) found that the Mach number distribu- 
tion peaks in the range 1 ^ Al < 3. Since diff'usive shock acceler- 
ation of CRs depends sensitively on the Mach number, this implies 
a more efficient CR injection in the simulations by Miniati et al. 
(2001a). (2) The grid-based cosmological simulations have been 
performed in a cosmological box of side-length 50 Mpc with 
a spatial resolution of 200 kpc, assuming an Einstein-de Sitter 
cosmological model (Miniati et al. 2001a). The lack of resolution 
in the observationally accessible, dense central regions of clusters 
in the grid-based approach underestimates CR cooling processes 
such as Coulomb and hadronic losses. Secondly, these simulations 
are unable to resolve the adiabatic compression of a composite of 
CRs and thermal gas, an effect that disfavours the CR pressure rela- 
tive to the thermal pressure. To summarise, their modest resolution 
in non-radiative simulations anticipates some of the results that we 
obtained using high-resolution simulations with radiative hydrody- 
namics and star formation, however for diff'erent reasons. 



4.2 Limitations and future work 

An accurate description of CRs should follow the evolution of the 
spectral energy distribution of CRs as a function of time and space, 
and keep track of their dynamical, non-linear coupling with the hy- 
drodynamics. We made several simplifying assumptions to enable 
the task of following CR physics self-consistently in cosmologi- 
cal simulations of smoothed particle hydrodynamics (SPH). In the 
following, we outline the possibly most severe limitations of our 
approach for computing the non-thermal emission processes that 
will be addressed in future work (cf. EnBlin et al. 2007, for a more 
complete list of the assumptions of our CR formalism). 

(i) We assumed a simple scaling of the magnetic energy density 
with the thermal energy density that allows us to eff'ectively scan 
the observationally allowed parameter space for the magnetic field. 
Note that current SPH implementations that are capable of follow- 
ing magneto-hydrodynamics (MHD) are presently still fraught with 
numerical and physical difficulties, in particular when following 
dissipative gas physics (Dolag et al. 1999, 2005; Price & Monaghan 
2004, 2005). While the inverse Compton and pion decay emission 
are mostly independent of the magnetic field, our synchrotron maps 
might be modified when the magneto-hydrodynamics is properly 
accounted for. 

(ii) We neglected the population of re- accelerated electrons 
throughout this work: strong merger shocks and shear motions at 
the cluster periphery might inject hydrodynamic turbulence that 
cascades to smaller scales, feeds the MHD turbulence and even- 
tually might be able to re-accelerate an aged CR electron popula- 
tion. Due to non-locality and intermittency of turbulence, this could 
partly smooth the very inhomogeneous primary emission compo- 
nent predominantly in the virial regions of clusters where simula- 
tions indicate a higher energy density in random motions. How- 
ever, to study these eff'ects, high-resolution AMR simulations are 



required that refine not only on the mass but also on some tracer 
for turbulence such as the dimensionless vorticity parameter. 

(iii) In our model, the emphasis is given to the dynamical im- 
pact of CRs on hydrodynamics, and not on an accurate spectral 
representation of the CRs. The pion decay emission is almost in- 
dependent on the spectral CR properties. However, the secondary 
CR component starts to be aff'ected by this simplification since the 
dimensionless CR momentum q ^ IGye^e/'^p - 100, that gives 
rise to synchrotron/IC emitting electrons with a Lorentz factor of 
7e ^ 10"^ is already quite high. Improving the spectral description 
of CR physics will not only allow us to study the spectral variations 
of the CR proton component but also enable reliable predictions 
for the TeV y-ray emission. This is of great interest for imaging air 
Cerenkov telescopes. 

(iv) We neglected microscopic CR diffiision in our simulations. 
The diffiisivity can be rewritten into a macroscopic advection term 
that we fully resolve in our Lagrangian SPH simulations by con- 
struction and a microscopic diff'usivity. The advection term domi- 
nates over microscopic term, in particular for merging clusters that 
are relevant for radio halos as the following estimate for the diffii- 
sivities shows: /Cadv = 100 kpc 1000 km/s = 10^^-^ cm^/s » Km = 
10^9 cmVs. 

(v) Our model of the diff'usive shock acceleration mechanism 
assumes a featureless power-law for both, the proton and the elec- 
tron acceleration, that is injected from the thermal distribution. The 
complete theoretical understanding of this mechanism is currently 
an active research topic that includes thermalization processes of 
the time evolution of the kinetic distribution of particles (Wolfe & 
Melia 2006) as well as non-linear eff'ects and magnetic field ampli- 
fication (Vladimirov et al. 2006). Phenomenologically however, we 
believe that there are strong indications for the diff'usive shock ac- 
celeration mechanism to be at work which come from observations 
of supernova remnants over a wide range of wavelengths from the 
radio. X-rays into the TeV y-rays (e.g., Ellison 2000; Hughes et al. 
2000; Ellison & et al. 2005; Warren et al. 2005; Aharonian et al. 
2004, 2006) as well as the bow shock of the Earth (Ellison et al. 
1990; Shimada et al. 1999). Future work will be dedicated on im- 
proving our model to incorporate more elaborate plasma physical 
models and to study the uncertainty of our results with respect to 
the saturated value of our CR acceleration efficiency (e.g., Kang & 
Jones 2007; Edmon et al. 2007). Varying the physics in our sim- 
ulations (non-radiative versus radiative) results in a very diff'erent 
Mach number distribution and changes the injection efficiency dra- 
matically (Paper I). However, the resulting non-thermal emission is 
almost independent of the simulated physics. For this reason, we 
are confident that our model produces reliable results due to the 
self-regulated nature of CR proton feedback. 

(vi) In this work, we did not account for feedback processes by 
AGN despite their importance for understanding the nature of the 
very X-ray luminous cool cores found in many clusters of galax- 
ies. In particular we neglected an additional CR population that 
diffiises out of AGN-inflated bubbles and postpone their study to 
future work (Sijacki et al. 2007). 

(vii) We furthermore neglected a hypothetical source of sec- 
ondary electrons that are produced "in-situ" in dark matter (DM) 
neutralino annihilations (Colafrancesco et al. 2006). However, we 
see two main challenges associated with the DM model for giant 
radio halos. (1) Only a sub-class of clusters that seems to be asso- 
ciated with past or present merging events exhibit radio halos. This 
fact makes it hard to believe that the dark matter with its nearly 
universal density profile in halos (Navarro et al. 1996, 1997) could 
be responsible for such an infrequent cluster property, especially 
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given the stringent upper limit on the diffuse radio emission in some 
massive clusters (Brunetti et al. 2007). (2) In order to explain the 
extended radio halo emission, the DM model needs to invoke a pro- 
file for the magnetic energy density, that increases by two orders of 
magnitude beyond the thermal core causing the plasma beta param- 
eter to decrease by a factor of 220 to a value of 15. This behaviour 
is not only in contrast to the magnetic profile predicted by numeri- 
cal MHD simulations of galaxy clusters Dolag et al. (2001) but also 
in contrast to turbulent dynamo models for the growth of the mag- 
netic fields that will saturate on a level which is determined by the 
strength of the magnetic back-reaction (e.g., Subramanian 2003; 
Schekochihin & Cowley 2006) and is typically a fraction of the tur- 
bulent energy density. Observationally, it is clear that the rotation 
measure values towards radio galaxies show a larger dispersion for 
smaller projected separation from the cluster centre in a large rota- 
tion measure sample (Clarke et al. 2001). This is mainly due to the 
distribution of these sources within the cluster atmosphere causing 
some foreground objects to experiencing a small Faraday rotation 
depth which is however no argument in favour of a decreasing mag- 
netic profile towards the cluster centre. 



5 CONCLUSIONS 

We find that the cosmic ray (CR) proton pressure traces the time 
integrated non-equilibrium activities of clusters and is only mod- 
ulated by the recent dynamical activities. In contrast, the pressure 
of primary shock accelerated CR electrons resembles the current 
dynamical, non-equilibrium activity of forming structure and re- 
sults in an inhomogeneous and aspherical spatial distribution. This 
is due to the diff'erent cooling time scales of CR electrons and pro- 
tons and reflects their large mass diff'erence. Hence, the radio syn- 
chrotron and inverse Compton emission of primary electrons pro- 
vides a snapshot of violent non-equilibrium processes that are re- 
sponsible for dissipating gravitational energy associated with struc- 
ture formation such as merger shock waves. Signatures of this emis- 
sion component are irregular morphologies, large spectral varia- 
tions, and a high degree of synchrotron polarisation. On the other 
hand, non-thermal emission processes of pions and secondary CR 
electrons produced in hadronic CR proton interactions trace the 
comparably smooth CR proton distribution centred on the cluster 
that the CR protons accumulate over the Hubble time. 



5.1 Radio synchrotron emission 

Unified scheme: we propose a unified scheme for the generation 
of giant radio halos, radio mini-halos, and radio relics that natu- 
rally arises from our simulated radio synchrotron maps and emis- 
sion profiles. It is schematically shown in a cartoon in Fig. 17. It 
predicts that the diff'use radio emission from a cluster varies with 
its dynamical stage as follows: 

(i) Once a cluster relaxes and develops cool core, a radio mini- 
halo develops due to synchrotron emission of hadronically pro- 
duced CR electrons. Adiabatic compression of magnetic fields dur- 
ing the formation of the cool core should confine the observable ra- 
dio synchrotron emission to the cooling region of the cluster. Since 
the cooling gas accretes onto the central black hole, this triggers the 
radio mode feedback of the AGN. Radio emission from the jet typi- 
cally outshines the diff'use mini-halo which implies a high dynamic 
flux range. This leads to a negative selection efl'ect that disfavours 
the detection of radio mini-halos or makes it very challenging. 



Unified scheme for diffuse cluster radio emission 
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Figure 17. A cartoon showing our unified scheme for the diffuse cluster 
emission. Hadronic cosmic ray proton interactions with ambient protons 
of the thermal ICM are thought to be responsible for radio mini-halos in 
cooling core clusters and the central parts of giant cluster radio halos in 
merging clusters. This emission mechanism produces a regular cluster-wide 
morphology resembling the thermal bremsstrahlung emission. In contrast, 
the radio synchrotron emission of shock accelerated cosmic ray electrons 
(through the Fermi 1 mechanism) is believed to be responsible for giant 
radio relics that have an inhomogeneous morphology and are primarily lo- 
cated at the cluster periphery. 



(ii) If a cluster experiences a major merger, two leading shock 
waves are produced at first core passage that move outwards and 
become stronger as they break at the shallow peripheral cluster po- 
tential. Relativistic electrons are efficiently accelerated by means 
of diff'usive shock acceleration and magnetic fields are amplified 
by shock compression and MHD turbulence at these shocks. Due 
to the short cooling times of r ~ 10^ yrs, the synchrotron radi- 
ating electrons are confined to a narrow emission volume around 
the shock wave. In combination with the preferred magnetic field 
direction in the shock surface, this implies a high degree of syn- 
chrotron polarisation. The observer will typically observe one or 
two large-scale radio relics, depending on the merger character- 
istics such as the relative cluster masses, concentrations, and gas 
fractions, the merger geometry with respect to the line-of- sight, as 
well as the time dependent merger stage. Our cosmological simu- 
lations supports the picture put forward in isolated cluster merging 
simulations (Roettiger et al. 1996). 

(iii) Simultaneously, virialisation of the gravitational energy, 
that is associated with the merger, generates a morphologically 
complex network of shock waves. The lower sound speed in the 
cluster outskirts imply stronger shocks that accelerate a spatially 
irregular distribution of CR electrons in these regions. The injected 
MHD turbulence amplifies magnetic fields through strong shear 
motions and turbulent dynamo processes. The induced radio syn- 
chrotron emission traces these non-equilibrium processes similarly 
as the water 'gischt' traces breaking non-linear waves. A giant ra- 
dio halo develops due to (1) boost of the hadronically generated 
radio emission in the centre and a transition to the (2) irregular 
radio 'gischt' emission in the cluster outskirts that represents ra- 
dio synchrotron radiation emitted from primary, shock- accelerated 
electrons. 

Predictions: The observational consequences of our unified 
scheme can be summarised as follows: 
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(i) Clusters undergoing major mergers are expected to have a gi- 
ant radio halo with an extended radio synchrotron emission region 
(^haio - 0.5/?2oo) while relaxed cool core clusters should host a 
smaller radio-mini halo. 

(ii) Our simulated radio luminosities reproduce observed lumi- 
nosities of halos/relics for magnetic fields derived from Faraday 
rotation measurements (cf. Paper III). 

(iii) The regular morphology of the central parts of giant radio 
halos is a consequence of the dominant contribution of hadronically 
produced electrons. 

(iv) The morphology and the radio spectral index in the radio 
halo periphery is predicted to show large variations due to the dom- 
inant contribution from primary CR electrons generated by shock 
waves. Superposing in projection many causally disconnected syn- 
chrotron emitting shock regions with diff'erent shock strength and 
thus electron spectral indices leads to spectral variations. 

(v) The amount of the primary radio emission depends critically 
on the characteristics of the merger. We thus expect a large scatter 
in the scaling relation of the radio halo luminosities with cluster 
mass as well as in the pixel-to-pixel correlation of the thermal X- 
ray brightness with radio surface brightness. 

(vi) The central radio emission should be Faraday depolarised 
assuming statistically isotropic distribution of magnetic field while 
the external emission regions are expected to have a small degree 
of polarisation. As a word of caution, in order to detect this polar- 
isation one might be forced to go out to large impact parameters 
where the resulting synchrotron surface brightness is small and the 
emission is dominated by very few contributing emission regions 
along the line-of- sight. 

These predictions from our cosmological high-resolution sim- 
ulations successfully reproduce characteristics of observed radio 
relics, giant radio halos, as well as radio mini halos (Feretti et al. 
2004; Cassano et al. 2006). In our approach, we choose the mag- 
netic energy density to scale with the thermal energy. The formation 
of a cool core is expected to compress the magnetic field adiabati- 
cally and should be responsible for the peaked central radio mini- 
halo emission profile. This eff'ect should reinforce our observed dif- 
ference in emission size between giant radio halos and mini-halos. 

The observed correlation between radio halos and merging 
clusters implies a departure of these systems from hydrostatic equi- 
librium and leads to a complicated non-spherical morphology. The 
resulting X-ray mass estimates are subject to large uncertainties 
which makes the analysis and theoretical model building based 
on azimuthally averaged quantities questionable if not impossible 
since it causes loss of information and might yield biased results. 
For this reason, we analyse pixel-to-pixel correlations of the ther- 
mal X-ray surface brightness with non-thermal cluster emission 
processes. We find that the hadronically induced non-thermal emis- 
sion is tightly correlated with the thermal bremsstrahlung emission 
with the slope depending on the realisation of the magnetic field. In 
contrast, the correlation is much weaker and the scatter is increased 
in the case of primary non-thermal emission where structures in 
the correlation space density correspond to individual structure for- 
mation shock waves. This implies that in general, simulations will 
only be able to provide quantitative predictions for the statistical 
behaviour rather than deterministic predictions for the correlations. 
Our new radio halo model matches qualitatively the observed tight 
correlation at high surface brightness which broadens and flattens 
towards dimmer brightness (cf. Govoni et al. 2001). However, our 
correlation is slightly steeper than observed ones. Taking into ac- 



count the uncertainty of the magnetic field model we conclude that 
observed pixel-to-pixel correlations support our model. 
Future: What can we learn from a future, large sample of clusters 
that show diffiise radio emission? 

(i) Radio relics and giant radio halos occur in dynamically merg- 
ing clusters and indicate a departure of these systems from hydro- 
static equilibrium and spherical symmetry. This has to be taken into 
account in the derivation of the cluster mass. 

(ii) The orthogonal information about the dynamical cluster ac- 
tivity, that in general can not be obtained from the thermal cluster 
observables such as X-ray emission and Sunyaev-Zel'dovich efi'ect, 
will help us in constructing a 'gold cluster sample' for cosmology. 

(iii) The property of the spatially confined radio relic emission 
from shock accelerated electrons might be employed to solve the 
inversion problem of reconstructing the course of a merger event 
given the thermal and radio synchrotron observables. 

(iv) Combining high-resolution X-ray, Sunyaev-Zel'dovich, and 
radio observations will allow us to probe fundamental plasma 
physics: diffusive shock acceleration, large scale magnetic fields, 
and turbulence. 

5.2 Inverse Compton and y-ray emission 

In principle, inverse Compton (IC) emission and high-energy y-ray 
emission from decaying pions, produced in hadronic CR interac- 
tions, is the cleanest way of probing current structure formation 
shock waves as well as time integrated non-equilibrium cluster ac- 
tivity. This is because these non-thermal emission components are 
not weighted with the magnetic energy density as it is the case for 
synchrotron emission. Our main findings can be summarised as fol- 
lows. 

(i) We identify two main regions for the generation of non- 
thermal emission in clusters of galaxies: the core that is also emit- 
ting thermal X-rays and the virial regions where the accretion 
shocks reside and merging shock waves break at the shallower clus- 
ter potential. 

(ii) In the cluster core regions, the emission for energies Ey > 
100 MeV is dominated by pion decay y-rays. At lower energies, 
the IC emission from secondary CR electrons dominate the emis- 
sion. Only in merging clusters, the situation may be reversed for 
the outer cluster regions where the primary IC emission can attain 
a similar flux level as the pion decay emission and even exceed the 
secondary IC emission at lower energies. 

(iii) While the total high-energy y-ray emission is always domi- 
nated by the pion decay component irrespective of the clusters dy- 
namical state, the total hard X-ray IC emission can be dominated 
by either primary or secondary emission components, depending 
whether a major merger takes place that boosts the primary IC 
emission. 

(iv) A corollary of this is that the high-energy y-ray emission 
can be reliably predicted for massive clusters using a scaling re- 
lation of non-thermal emission and the cluster mass. In contrast, 
the hard X-ray emission of even massive clusters is subject to large 
flux variations that depend sensitively on the dynamical state of the 
cluster. 

(v) Due to larger variation of merging histories and the smaller 
gravitational potential in less massive systems, their CR pressure 
and the associated y-ray emission level is subject to larger modu- 
lation and reflects more sensitively the current merging activity of 
the cluster than it is the case in large systems. 

(vi) The morphology of the pion decay as well as the secondary 
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IC component resemble the thermal X-ray emission albeit they de- 
crease less steeply with growing radius and extend further out. This 
is due to the increasing CR number fraction /cr = ^cR/^th with in- 
creasing radius and reflects the more efficient CR acceleration at 
stronger shocks in the cluster periphery. The morphology of the 
primary IC emission is irregularly shaped and traces current non- 
equilibrium phenomena such as merger or accretion shock waves. 

(vii) Possibly most surprising, we find that the dominant emis- 
sion component at the centre (primary or secondary IC for Ey > 
10 keV and pion decay y-rays for Ey > 100 MeV) depends only 
weakly on whether radiative or non-radiative gas physics is simu- 
lated provided we consider in both cases only CRs from structure 
formation shocks. This is mainly due to self-regulating effects of 
the CR pressure. 

(viii) Measuring the hard X-ray and y-ray emission will have a 
huge astrophysical impact and teach us about: the CR pressure con- 
tribution to the intra-cluster medium, the generating mechanisms 
of radio halos (such that we can use them in addition to thermal 
observables to characterise clusters), the contribution of the pion 
decay emission as well as the primary and secondary IC radiation 
to the y-ray background. 

(ix) Detecting the non-thermal spectrum ranging from X-rays 
to y-rays will enable us to probe fundamental plasma physics on 
large cluster scales such as inferring the energy conversion effi- 
ciency of diff'usive shock acceleration of protons and electrons as 
well as probing the large scale magnetic fields. 
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APPENDIX A: RELATIVISTIC ELECTRON 
POPULATIONS 

Al Definitions 

Throughout the paper we use the following definitions for the dif- 
ferential source function q{r, E), the emissivity j(r, E) and the vol- 
ume integrated quantities, respectively: 



q{r,E) = 
Q(E) = 



dtdVdE 

J dVq{r, E), 



j(r,E) = Eq(r,E), (Al) 
J(E) = EQ(E), (A2) 



where N denotes the integrated number of particles. From the 
source function the integrated number density production rate 
of particles A(r), the number of particles produced per unit time 
interval within a certain volume, X, and the particle flux 9^ can 
be derived. The definitions of the energy weighted quantities are 
denoted on the right hand side, respectively. 



£ 

r 



J dEq(r,E), 
J dVA(r), 



A(r) 
L 
F 



J dEEq(r,E), (A3) 
J dVA(r), (A4) 

4^- 



A2 Timescales 

This section presents general considerations for deriving the char- 
acteristic electron momentum scales of the distribution function. To 
this end, we compare the energy loss timescales tioss = -E/Eio^^ 
to the acceleration timescales as a function of particle kinetic en- 
ergy for the most important processes in this context. The de- 
tailed calculations for the equipartition distribution of CR electrons 
can be found in the next two sections. The particle kinetic energy 
EKiUqC^) = y- 1 = (1 -0^)~^l^ - 1 is defined in terms of the Lorentz 
factor 7 and the velocity of the electron, ^c. Energy loss processes 
due to Coulomb interactions (Gould 1972b), inverse Compton (IC) 
and synchrotron emission (Rybicki & Lightman 1979) are defined 
as follows: 



3 (Tt MqC^ 



In 



y- 



Ay 



-E, 



IC,synch — 



(A6) 
(A7) 

Here ctt - is the Thomson cross section, = 6^/(171^0^) 

the classical electron radius, co^\ = ^j4ne^nQ/mQ is the plasma 
frequency, and is the number density of free electrons. IC 
losses of electrons depend on the energy density of the cosmic mi- 
crowave background (CMB) and the starlight photon field, = 
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Figure Al. Energy gain and loss timescales as a function of the kinetic en- 
ergy of electrons for typical conditions of the ICM. The three solid lines 
from the bottom to the top denote the total loss timescale for electrons, 
the timescale due to hadronic injection of secondary electrons (pp), and the 
combined inverse Compton (IC)/synchrotron cooling timescale. The dotted 
line shows the Coulomb timescale, the long dashed one the synchrotron 
timescale, and the dashed one the Hubble time. Note that the relative nor- 
malisation of the hadronic injection timescale compared to the total loss 
timescale is subject to different density dependencies and the shock accel- 
eration timescale depends crucially on the properties of structure formation 
shocks. 



^CMB + ^stars. whcrc wc ucglcct the latter one for simplicity and ex- 
press £cMB = ^cmb/(^^) equivalent field strength Bqmb = 
3.24 (1 + z)^\xG. Synchrotron losses of an isotropic pitch angle dis- 
tribution of electrons depend on the energy density of the local 
magnetic field, sb = {B^)/(S7i), where B = ^J{B^) is the rms of 
the magnetic vector field B. Comparing these two loss processes, 
we obtain a synchrotron dominated cooling regime for B > Bqmb 
and an IC dominated regime in the weak field limit. The timescale 
for CR electron injection by means of hadronic interactions of CR 
protons with ambient protons of the thermal plasma is given by: 



-£'pp = Eccr^^n^, 



(A8) 



where o-pp = 32 mbarn is the average inelastic cross section for 
proton-proton interactions, = + Ariue = p/^p is the number 
density of target nucleons. 

Comparing the diff'erent energy gain and loss rates of rela- 
tivistic electrons yields characteristic momentum scales that are re- 
sponsible for spectral breaks or cutoff's in the CR electron distri- 
bution function. Conveniently, we denote these characteristic mo- 
mentum scales with the dimensionless electron momentum q = 
jSy = Pe(mec)~^ and label these with the competing processes 
considered. Equating the timescale for Coulomb interactions and 
IC/synchrotron losses yields 



^Coul,IC/synch 




(^!!^]. 0.216 



,1/2 



(A9) 
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Figure A2. Density dependence of characteristic momentum scales q = 
Pq (mec)~^ for a hadronically injected secondary CR electron population. 
Comparing the different energy gain and loss rates of competing processes 
such as hadronic injection of secondary electrons (pp), Coulomb cool- 
ing, and IC/synchrotron cooling yields characteristic momentum scales that 
are responsible for spectral breaks or cutoffs in the CR electron distribu- 
tion function. The solid red line denotes the effective spectral cutoff: for 
«e 10~^ cm~^, it denotes the lower cutoff of the power-law distribution 
function. At higher densities, it denotes the lower cutoff ^couijc and the 
spectral break (^ppjc) of the broken power-law spectrum. The dotted lines 
include synchrotron losses and assume a scaling of the magnetic field of 
°^ ^th with a saturation value for B at -Smax = 10 



In the last step and the following numerical examples, we assume 
the IC cooling regime for simplicity. This energy scale provides a 
bottleneck through which high-energy electrons have to pass when 
they age and one expects a characteristic pile-up at this energy scale 
in their distribution function for an integration over the energy in- 
jection history of different CR electron populations (Sarazin 1999). 

The timescale for diffusive shock acceleration of electrons at 
cosmological shock waves is much shorter than cosmologically rel- 
evant timescales, Tshock ^ 1 Gyr. This implies that Tshock inter- 
sects the total loss timescale in the low-energy Coulomb regime, 
Tcouh as well as in the high-energy IC/synchrotron regime. Tic/synch 
(cf. Fig. Al). We can thus identify two characteristic momenta of 
the primary population of electrons that are obtained by equating 
Tshock with the Coulomb and the IC/synchrotron timescale. 



^inj, IC/synch 
^inj,Coul 



3 m^c 



4 CTxTinj {Sb + fiph) 



- CTxC^eTinj 



in, 0.216 



'pi 



(AlO) 
(All) 



At high energies, we expect to have an IC/synchrotron cooled elec- 
tron spectrum that joins at lower energies into the shock injection 
spectrum of CR electrons that have had no time to cool radiatively 
yet. The low momentum cutoff of the CR electron distribution func- 
tion is determined by Coulomb losses. 

Characteristic momenta of the secondary population of elec- 
trons are obtained by equating the energy injection rate through 
hadronic proton interactions with the energy loss rates. 



^ 70, 



3,IC/synch 



4 me 

3 d-pp n^ me 
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(A13) 
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2 O-pp «N 

^ 1300 - In 
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(A14) 



^threshold reflects the threshold momentum for the inelastic proton- 
proton reaction. The shortest equipartition timescale at these 
characteristic momenta dominates the resulting electron equilib- 
rium distribution. In the case of a double-valued solution for 
q, i.e. when we obtain two equipartition cutoffs with a similar 
timescale, we choose the larger one which is in equilibrium with 
the IC/synchrotron cooling. All these momentum scales have dif- 
ferent density dependencies which are visualised in Fig. A2. At 
typical densities of the ICM for n^ < 10~^ cm~^, the two momen- 
tum scales ^ppjc/synch and ^pp,coui have an associated timescale that 
is much longer than the timescale at ^couijc/synch as can be readily 
inferred from Fig. Al. This implies that the equilibrium distribu- 
tion function of secondary CR electrons has a low-energy cutoff at 

^CoulJC/synch- 

The situation is reversed for the dense ICM or the inter- stellar 
medium with electron densities n^ > 10~^ cm~^ (assuming the IC 
dominated cooling regime), and the momentum scale ^couijc/synch 
drops out of the problem due to its long timescale. The equilibrium 
distribution function of secondary CR electrons develops a break at 
^pp,ic/synch abovc which the secondary electron injection is in equi- 
librium with IC/synchrotron cooling and below which the electron 
injection spectrum remains unchanged (similar to the case of the 
primary CR electron population). The lower cutoff of the distribu- 
tion function is given by ^pp,coub provided it exceeds the threshold 
^threshold for the hadronic reaction and provided the CR proton dis- 
tribution extends down to these low energies. The effective spec- 
tral cutoff ^eff of the electron distribution function is visualised in 
Fig. A2 as a solid red line. For n^ < 10~^ cm~^, we have a simple 
power-law with a lower cutoff given by q^s. At higher densities, 
the lower solid red line denotes again the cutoff of the distribu- 
tion function, while the upper red line (that coincides with ^ppjc) 
denotes the spectral break above which the hadronic injection and 
IC/synchrotron cooling established a steady state spectrum. 

So far, we only considered regime of weak magnetic fields 
where IC cooling dominates. The complete picture including syn- 
chrotron cooling is however only slightly changed due to the fol- 
lowing line of arguments. Assuming a simple scaling model for the 
magnetic energy density, Sb oc n^, in Eqns. (A9) and (A 13) will 
cause these momentum scales to become independent of density, 
as modelled in Fig. A2 with dotted lines. Eventually, the enhance- 
ment of the magnetic field strength through turbulent dynamo pro- 
cesses will saturate on a level which is determined by the strength 
of the magnetic back-reaction (e.g., Subramanian 2003). For il- 
lustrative purposes in Fig. A2, we model such a saturation effect 
with a simple model, sb = sbq (i - exp(«e/«eo))^ where we chose 
^max = 10 \iG and n^^ = 10~^ cm~^. This causes the transition to 
the broken power-law CR electron spectrum to occur at the higher 
critical electron density 



ne.cnt = 10 cm 

^CMB 



(A15) 



where sb^^^ - 5^^/(87r) and we assumed a plasma of primor- 
dial composition with a hydrogen mass fraction of Xh = 0.76 
and full ionisation. These modifications due to synchrotron losses 
cause a curvature of the straight lines in the log-log representation 
in Fig. A2 of the original power-law dependencies on n^. As before, 
shown in red is the spectral cutoff and spectral break, respectively. 
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above which the hadronic injection and IC/synchrotron cooling es- 
tablished a steady state spectrum. 



A3 Shock-accelerated electron population 

A3.1 Injection spectrum 

In this section, we discuss electron acceleration processes at shock 
waves due to gas accretion and galaxy mergers, using the frame- 
work of diffusive shock acceleration using the thermal leakage 
model originally proposed by Ellison et al. (1981). Our description 
follows the approach of EnBlin et al. (2007) for the acceleration of 
CR protons. The shock surface separates two regions: the upstream 
regime defines the region in front of the shock which is causally 
unconnected for super sonic shock waves, whereas the wake of the 
shock wave defines the downstream regime. The shock front itself 
is the region in which the mean plasma velocity changes rapidly 
on a scale of the order of the plasma skin depth. In the rest frame 
of the shock, particles are impinging onto the shock surface at a 
rate per unit area of P2V2 = Pi^i- Here vi and V2 give the plasma 
velocities (relative to the shock's rest frame) in the upstream and 
downstream regimes of the shock, respectively. The corresponding 
mass densities are denoted by pi and p2. 

We assume that after passing though the shock front most of 
the electron gas thermalizes to a Maxwell-Boltzmann distribution 
with characteristic post-shock temperature T2 and the dimension- 
less electron momentum p = Pq {m^cY^'. 



Mp) = 4nn, 



m^c^ 
2nkT2 



3/2 



p^ exp 



2 2 
m^c p 

2kT2 



(A16) 



where the number density of electrons of the thermal distribution 
in the downstream regime, n^ = ^2, as well as T2 can be in- 
ferred by means of the mass, momentum, and energy conservation 
laws at the shock surface for a gas composed of relativistic parti- 
cles and thermal constituents. Note that we use an eff'ective one- 
dimensional distribution function f(p) = 4np^f^^\p). In our sim- 
ulations, we follow the spatial and temporal evolution of the hy- 
drodynamic quantities such as temperature and density (although 
for brevity we suppress this in our notation). For cosmological ap- 
plications, we have to consider the primordial composition of the 
cosmological fluid, i.e. the ionised electron number density is given 
by We = XnXep/mp, where Xu = 0.76 is the primordial hydrogen 
mass fraction, and is the ratio of electron and hydrogen number 
densities which we dynamically track in our radiative simulations. 
Assuming that a fraction of the thermalized particles experience 
stochastic shock acceleration by diff'using back and forth over the 
shock front, the test particle theory of diffiisive shock acceleration 
predicts a resulting CR power-law distribution in momentum space. 
Within our model, this CR injection mechanism can be treated as 
an instantaneous process. 

For a particle in the downstream region of the shock to re- 
turn upstream it is necessary to meet two requirements. The par- 
ticle's eff'ective velocity component parallel to the shock normal 
has to be larger than the velocity of the shock wave, and secondly, 
its energy has to be large enough to escape the "trapping" process 
by Alfven waves that are generated in the downstream turbulence 
(Malkov & Volk 1995; Malkov & Volk 1998). Thus, only parti- 
cles of the high-energy tail of the distribution are able to return 
to the upstream shock regime in order to become accelerated. The 
complicated detailed physical processes of the specific underlying 
acceleration mechanism are conveniently compressed into a few 
parameters (Jones & Kang 1993; Berezhko et al. 1994; Kang & 



Jones 1995), one of which defines the momentum threshold for the 
particles of the thermal distribution to be accelerated. 



^inj ~ ^mjPe ~ -^in 



(A17) 



(A19) 



In the linear regime of CR electron acceleration, the thermal 
distribution joins in a smooth manner into the resulting CR elec- 
tron power-law distribution at ginj so that xinj represents the only 
parameter in our simplified diffusive shock acceleration model, 

Mp) = /e(^inj) (^) 0(p - ^inj). (A18) 

The slope of the injected CR electron spectrum is given by 

r + 2 Pi y\ 

<3^ini = 7 , where r = — = — 

r - 1 P\ vi 

denotes the shock compression ratio (Bell 1978a,b; Drury 1983a). 
In combination with the slope ofinj, the value of xinj regulates the 
amount of kinetic energy which is transferred to the CR electrons. 
Theoretical and observational studies of shock acceleration of CR 
protons at galactic supernova remnants suggest a range of xinj ^3.3 
to 3.6, implying a particle injection efficiency of T/un - 10""^ to 10~^ 
(Drury et al. 1989; Jones & Kang 1993; Berezhko et al. 1994; Kang 
& Jones 1995; Malkov & Volk 1995; Ellison 2000; Hughes et al. 
2000; Warren et al. 2005). In the linear regime, the number density 
of injected electrons is given by 



^n^i, 



Jo 



dpflinip) = /e(^inj) 



^inj 



1 



(A20) 



This enables us to infer the particle injection efficiency which is a 
measure of the fraction of downstream thermal gas particles which 
experience diffusive shock acceleration. 



^lin 



n. 



^/n Q^inj - 1 



(A21) 



The particle injection efficiency is independent of the downstream 
post- shock temperature T2. These considerations allow us to infer 
the injected electron energy density in the linear regime: 



Aeiin = ^lin^r^(a^inj,^inj)'^e(7'2). 



(A22) 



The average kinetic energy of ^^e^^Q^inj, ^inj) of an injection 
power-law spectrum with CR electron spectral index ainj and lower 
momentum cutoff ^inj is given by 



— = «r„j r dp^j(p)£(p) 

^inj Jo 

q""'^ ^ la-2 3-a\ r ^ 



XA23) 



where E{p) = ( ^/l + - 1) me c^ is the kinetic energy of an 
electron with momentum p, the electron distribution function /inj 
is given by Eqn. (A27), and we used the abbreviations a = 
CKinj and q = qiny Bx{a,b) denotes the incomplete beta-function 
(Abramowitz & Stegun 1965), assuming a > 2. In our descrip- 
tion, the CR electron energy injection efficiency in the linear regime 
is defined to be the energy density ratio of freshly injected CR 
electrons to the total dissipated energy density in the downstream 
regime. 



^in 



A£iii 
ASdi! 



where Asdiss = s^i - s^xr^ - 



(A24) 



The dissipated energy density in the downstream regime, A^dis 
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is given by the difference of the thermal energy densities in the 
pre- and post-shock regimes, corrected for the contribution of the 
adiabatic part of the energy increase due to the compression of the 
gas over the shock. 

In order to obey energy conservation as well as the saturation 
effect for strong shocks, we propose the following modification of 
the electron injection efficiency at high values of the Mach number: 



^inj — 



l-exp|-|i^ 

bmax 



(A25) 



Keshet et al. (2003) suggest a value of ^max - 0-05 for the limit- 
ing case of the electron energy injection efficiency. One can then 
infer the injected CR electron energy density in terms of the energy 
injection efficiency of diffusive shock acceleration processes. 



^^inj — ^inj^^diss- 



(A26) 



Putting these considerations together, one arrives at the injec- 
tion spectrum for the electrons, 



fmi(p)dp = Cinj P "^'"J e(p - ^inj)d/7, 

A^e.lin ^Q,lin Eq (o^inj , ^inj ) 



6 = 



/etanj) = 



^e,lin 



^max ^diss^shock 



4 



(A27) 
(A28) 

(A29) 
(A30) 



where ^diss = ^diss, sph^^p/C^sph^h -^e) denotes the dissipated en- 
ergy per timestep and per electron and Tshock = fhh/v is the time 
it takes the particle to pass through the broadened shock front. The 
front has a characteristic length scale that is a multiple of the SPH 
smoothing length h (with fh = 2), and one may approximate v with 
the pre-shock velocity vi = AtiCgi. 



A3. 2 Equilibrium spectrum of shock accelerated electrons 

This section describes the steady- state approximation for the equi- 
librium CR electron spectrum. This is only justified if the dynami- 
cal and diffusive timescales are long compared to the shock injec- 
tion or IC/synchrotron timescale. This may well be the case in clus- 
ters of galaxies, however, probably not in our own Galaxy. More- 
over, this section neglects possible re- acceleration processes of CR 
electrons like continuous in- situ acceleration via resonant pitch an- 
gle scattering by compressible MHD modes. 

The steady- state CR electron spectrum at high energies p = 
Py > GeV/c is governed by the injection of shock-accelerated CR 
electrons, denoted by the source function ^e, and their cooling pro- 
cesses so that it can be described by the continuity equation 



^ [P{P)UP)] 



Sc(p) . 



For p(p) < 0, this equation is solved by 



Up) = 



's^{p'). 



(A31) 



(A32) 



For the energy range of interest, the cooling of the radio emit- 
ting CR electrons is dominated by synchrotron and inverse Comp- 
ton losses, /?ic,synch = £ic,synch/(meC^) whcrc £ic,synch is given by 
Eqn. (A7). The source function of the shock-accelerated CR elec- 
trons for the energy range of interest is given by 

/inj(p) 



s^{p) = 



(A33) 



In our formalism, we set Tinj = min(Tshock,THubbie) using Tshock of 
Eqn. (A29) due to the following line of arguments: the freshly 
accelerated relativistic electron population in post-shock regions 
cools and finally diminishes as a result of loss processes. In the in- 
teresting observational bands such as inverse Compton (IC) y-rays 
and radio synchrotron emission the electron population diminishes 
on such a short timescale that we could describe this by instanta- 
neous cooling. In this approximation, there is no steady-state elec- 
tron population and we would have to convert the energy from the 
electrons to IC and synchrotron radiation. However, we can intro- 
duce a virtual electron population that lives in the SPH broadened 
shock volume only which is defined to be the volume of energy dis- 
sipation. Within this volume that is comoving with the shock, we 
can indeed use the steady state solution for the distribution function 
of relativistic electrons and we assume no relativistic electrons in 
the post-shock volume where there is no energy dissipation. Thus, 
the cooled CR electron equilibrium spectrum can be derived from 
Eqn. (A32) yielding 



Up) dp 



Cep-"^dp, 



Ce = 



3 Cinj m,QC 



4 (Ofe - 2) (Tt Tinj (^5 + ^ph) 



(A34) 
(A35) 



Here, = orinj + 1 is the spectral index of the equilibrium electron 
spectrum. The normalisation scales linearly with the gas density 
Cq oc p which we evolve dynamically in our simulations and de- 
pends indirectly on ofinj and £^diss through the variable Cjnj. 

At high energies, we have the IC/synchrotron cooled power- 
law electron spectrum that joins at lower energies into the shock 
injection spectrum which has had no time to cool radiatively yet. 
The low-energy regime of the CR electron distribution function is 
determined by Coulomb losses. It turns out that the timescale asso- 
ciated with the momentum scale ^couijc/synch is always larger than 
the injection timescale Tinj such that the transition from the injection 
spectrum (A27) to the cooled equipartition spectrum (A34) takes 
place at the characteristic momentum 



^break,prim — 



3 MqC 



4 {a^ - 2) CTxTinj {Sb + Sph) ' 



(A36) 



The injection spectrum extends down to the lower cutoff ^niin,prim = 
max(^inj,^inj,coui) where ^inj and ^inj,coui are given by Eqns. (A17) 
and (All), respectively. 

The pressure of a CR electron or proton power-law popula- 
tion as e.g. described by Eqn. (A34), that is characterised by two 
momentum cutoffs p\ and p2 is given by 

mc^ r 

PcR = — dpf(p)f3p (A37) 



Cmc^ 



a -2 3 - a 



(A38) 



where := v/c = p/ ^Jl + p^ is the dimensionless velocity of 
the CR particle. The CR population can hydrodynamically be de- 
scribed by an isotropic pressure component as long as the CRs are 
coupled to the thermal gas by small scale chaotic magnetic fields. 
Note, that for 2 < (2 < 3 the kinetic energy density and pressure of 
the CR populations are well defined for the limit q ^ 0, although 
the total CR number density diverges. 



A4 Hadronically produced electron population 

Considering CR protons, which are at least in our Galaxy the dom- 
inant CR species, it is convenient to introduce the dimensionless 
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proton momentum = Pp/(mp c). We assume that the differential 
particle momentum spectrum per volume element can be approxi- 
mated by a single power-law: 



dN 
dp^dV 



= CpP^ ^ e(pp - ^p), 



(A39) 



where 0(x) denotes the Heaviside step function, Cp = Cp(x, t) de- 
notes the normalisation, = qp(x, t) is the lower cutoff of the dis- 
tribution function, and ofp is the CR spectral index that is taken 
to be constant in space and time for simplicity. In our simulation, 
we dynamically evolve the quantities Cp and according to the 
dominant gain and loss processes in the intra-cluster medium. The 
modelling of the cosmic ray physics includes adiabatic CR trans- 
port processes, injection by supernovae and cosmological structure 
formation shocks, as well as CR thermalization by Coulomb inter- 
action and catastrophic losses by hadronic interactions. As already 
laid out in the introduction, the hadronic reaction of CR protons 
with ambient thermal protons produces pions which decay into sec- 
ondary electrons, positrons, neutrinos and y-rays. 

There are two analytical models in the literature that describe 
the hadronic proton-proton reaction while assuming isospin sym- 
metry. Fermi (1950) proposed the fireball model which assumes a 
state of hot quark-gluon plasma in thermal equilibrium after the 
hadronic interaction that subsequently ablates pions with energy 
dependent multiplicities. Since this model is only valid in the high- 
energy limit for CR protons, we use the analytic formalism by 
Pfrommer & EnBlin (2004a) that parametrises important effects 
near the pion threshold and is based on an approximate description 
developed by Dermer (1986a,b), which combines isobaric (Stecker 
1970) and scaling models (Badhwar et al. 1977; Stephens & Bad- 
hwar 1981) of the hadronic reaction. 



A4. 1 Injection spectrum 

The pion production spectrum can be derived from general con- 
siderations including branching ratios and multiplicities of the 
hadronic reaction (Stecker 1971). The pion production spectrum 
describes the produced number of pions per unit time, volume and 
momentum intervals, dN/(dt dV dp„dpp), and reads in this context 

Sn(Pn, Pp) = Cn^^(pp)0-l^(Pp)6B(p„ - {Pn))e{p^ " Pp,th), (A40) 

where = + 4nHe = p/i^^^ the target density of nucleons in a 
fluid of primordial element composition, cr^^ the inelastic p-p cross 
section, {p„) the average momentum of a single produced pion, and 
Pp,th = 0-78 denotes the threshold momentum for pion production. 
For a differential CR proton distribution, the pion source function 
can be marginalised over the proton energy, yielding 

2 2 

Sn^(Pn) = :^qn(Pn) ^ 3 J ^Ppfp(Pp)Sn(P7r, Pp), (A41) 

where the CR proton population is given by Eqn. (A39). The scal- 
ing behaviour in the high-energy limit of Dermer' s model can be 
described by a constant pion multiplicity ^(pp) ^ ^ = 2 and the 
dependence of the mean pion momentum is given by (pn^) = 
mpPp/(2mj,±^). The weak energy dependencies of the pion multi- 
plicity and the inelastic cross section can be absorbed in a semi- 
analytical parametrisation of the cross section, a^^ia^) (for details 
see Pfrommer & EnBlin 2004a). The mean energy of the produced 
secondary electrons (tt- e- + 3y) in the laboratory frame is given 
by (Eq) = \{E„±). Employing the transformation law for distribu- 
tion functions and using the mean value of the electron momentum 



in the relativistic Umit allows us to approximate the electron source 
function by 



4me I 






m„± \ 


ym„± 








dp. 



4 

= - 16^-^PCO-ppWN 

where the effective cross section cTpp depends in our model on the 
spectral index of the CRp spectrum according to 



o-pp ^ 32 . (0.96 + e4-4-2-4o^p) mbarn. 



(A44) 



Thus, we can write down the injection spectrum for CR electrons 
resulting from hadronic reactions of CR protons with ambient gas 
protons. 



/inj,pp dp = Cinj,ppP dp 

4 



(A45) 

(\2-ae 
^ , (A46) 
mp/ 

and Tpp = min[(ccrppnN)"\THubbie], and = + 1. 



C'inj,pp — 3 ^ ^ ^pp ^pp '^N C'p 



A4.2 Equilibrium spectrum of secondary electrons 

The same line of arguments presented in Sec. A3. 2 allows us to 
derive the equilibrium distribution of secondary CR electrons above 
a GeV due to IC and synchrotron cooling. 



Mp) dp 



Ce = 



C,p-^^ dp 

16^~"^^o-ppnNCp me c^ 

(Ofe - 2) O-T {Sb 



,me c^ /^r"^ 
+ %h)\me/ 



(A47) 



(A48) 



where the effective CR-proton cross section cTpp is given by 
Eqn. (A44), and n^ = + 4«He - p/^p is the target density 
of nucleons in a fluid of primordial element composition. As dis- 
cussed in Sect. A2, the equilibrium spectrum of secondary electrons 
looks different depending on the ambient electron density relative 
to the critical electron density (A 15). At average ICM densities be- 
low We,crit» the equilibrium spectrum is given by the IC/synchrotron 
cooled spectrum (A47) with the lower cutoff 



■ inax(^Coul,IC/synch, ^threshold » ^pp)» 



(A49) 



where ^pp = q^mp/ilGm^) is the lower cutoff of the injected elec- 
tron population that is inherited from the lower proton cutoff q^, 
while ^couijc/synch and ^threshold are given by Eqns. (A9) and (A12), 
respectively. 

In the high ICM densities/ISM-regime above ^ccrit, the equi- 
librium spectrum is given by the injection spectrum (A45) at low 
energies between ^niin,sec and 

^pp,IC/synch 



^break,sec 



(A50) 



Above ^break,sec, the cquipartition spectrum steepens and and joins 
continuously into the IC/synchrotron cooled spectrum (A47). 



APPENDIX B: RADIATIVE PROCESSES 

The non-thermal radio and hard X-ray emission is generated by 
CR electrons with energies E^ > GeV (cf. Eqns. (1) and (2)). For 
convenience, we rescale the cooled CR electron equilibrium spectra 
of Eqns. (A34) and (A47) to the energy scale of a GeV, 



UE,)dE, = 



Ce / E, \- 

GeV \GeV/ 



d^e, 



(Bl) 
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(2 \' 



2 \«'e-l 



(B2) 



and Ce is given by Eqn. (A35) respectively (A48), depending on 
the electron population. 



Bl Cluster magnetic fields 

In principle, cosmological structure formation calculations with 
SPH are capable of following magneto-hydrodynamics (Dolag 
et al. 1999, 2005; Price & Monaghan 2004, 2005), although this 
is presently still fraught with numerical and physical difficulties. 
Secondly, the origin of cluster magnetic fields is still an open ques- 
tion (Widrow 2002, and references therein). There are studies of 
the Faraday rotation measure (RM) as a function of cluster impact 
parameter using the position of a sample of radio lobes in diff'erent 
clusters (Clarke et al. 2001) which hints at a magnetic profile cen- 
tred on the cluster with |xG field strengths. Field reversals along the 
line-of-sight lead to cancellations in RM, since RM oc j n^B • 61. 
The unknown behaviour of the characteristic length scale of the 
magnetic field with cluster radius leaves us with some degree of 
freedom for the magnetic profile that is unconstrained by cur- 
rent observations. Assuming primordial origin, and amplification 
of magnetic fields in the process of structure formation would still 
require scanning the parameter space of the field strengths in the 
initial conditions (Dolag et al. 1999). 

Thus, we refrain from running self-consistent MHD simula- 
tions on top of the radiative gas and CR physics and postpone a 
detailed analysis of the influence of MHD on the radio emission to 
future work. We chose the following simple model for the magnetic 
energy density: 



gth 

\ ^th,0 



2aB 



(B3) 



where sb,o and or^ are free parameters in our model. Rather than 
applying a scaling with the gas density as non-radiative MHD sim- 
ulations by Dolag et al. (1999, 2001) suggest, we chose the en- 
ergy density of the thermal gas. This quantity is well behaved in 
the centres of clusters where current cosmological radiative sim- 
ulations, that do not include feedback from AGN, have an over- 
cooling problem which results in an overproduction of the amount 
of stars, enhanced central gas densities, too small central temper- 
atures, and too strong central entropy plateaus compared to X-ray 
observations. Theoretically, the growth of magnetic field strength is 
determined through turbulent dynamo processes that will saturate 
on a level which is determined by the strength of the magnetic back- 
reaction (e.g., Subramanian 2003; Schekochihin & Cowley 2006) 
and is typically a fraction of the turbulent energy density that itself 
should be related to the thermal energy density, thus motivating our 
model theoretically. 



B2 Synchrotron radiation 

The synchrotron emissivity jy at frequency v and per steradian of a 
CR electron population described by Eqn. (Bl), which is located in 
an isotropic distribution of magnetic fields (Eqn. (6.36) in Rybicki 
& Lightman 1979), is obtained after averaging over an isotropic 
distribution of electron pitch angles, yielding 



^%nch(«^e)Ce 



£b_ 



(«y+l)/2 



, Inmlc^ y ( v \ 



synch 



3271 rrieC^ a. + I 



r(^) 



(B5) 



(B6) 



where r(d) denotes the F-function (Abramowitz & Stegun 1965), 
ay = (ofe - 1)/2 = Qfinj/2, Ce is givcu by Eqn. (B2), and Be denotes a 
(frequency dependent) characteristic magnetic field strength which 
implies a characteristic magnetic energy density sb^ • Line-of-sight 
integration of the radio emissivity jy yields the surface brightness 
of the radio emission Sy. 

For later convenience, we calculate the radio luminosity per 
unit frequency interval of a steady state population of hadronically 
generated electrons (A47), 

/r ^ / ^ \(Q^y-i)/2 

dVjy = Ay I dVCp^N 



+ %h \ ^B^ 



SB 



^ Ay J dVCp^N, for sb » Sph, 



4nAE 



16 



2-ae 



CTpp me I mp \''' / meC^ V 
2)crTSB, \mj \GeV/ 



2 \0'e-l 



(B7) 



(B8) 



^synch 

where we introduced the abbreviation Ay with the dimensions 
[Ay] = erg cm^ s~^ Hz~^ and the volume integral extends over the 
entire cluster. In the last step of Eqn. (B7), we assumed typical ra- 
dio spectral indices of cluster radio halos of ofy ~ 1 such that the 
radio luminosity of the equilibrium distribution of CR electrons be- 
comes independent of the magnetic field in the synchrotron domi- 
nated regime for sb » Sph (cf. Fig. 3). 

B3 Inverse Compton radiation 

Inverse Compton (IC) scattering of cosmic microwave background 
(CMB) photons off" ultra-relativistic electrons with Lorentz factors 
of ye ~ 10"^ redistributes these photons into the hard X-ray regime 
according to Eqn. (2). The integrated IC source density Aic for an 
isotropic power law distribution of CR electrons as described by 
Eqn. (A34) or (A47), can be obtained by integrating the IC source 
function Sy(Ey) in Eqn. (43) of Pfrommer & EnBlin (2004a) (in 
the case of Thomson scattering) over an energy interval between 
observed photon energies Ei and E2 yielding 

dEic sic(Eic) 



Aic(Ei,E2) 



(B9) 



= ^o/lc(Q^e) 



/ me \ 


l-Q-e 


\i n 






\kTcMB) 



fic(ae) = 



and Aq = 



2^^+^ (ai + 4 a, + 11) 
(ore + 3)^ (ae + 5)(ae + I) 
I ae + 5\ I a. + 5] 



xF 



2 

16 TT^r^^ a 



(kTcMBr 



(BIO) 



(Bll) 



(B12) 



oc 6:cRe B"^^'y-^\ (B4) TyiE^E^) 



where ay = (a^-l)/! denotes the spectral index, - e^/(me c^) the 
classical electron radius, ^(d) the Riemann f-function (Abramowitz 
& Stegun 1965), and Q is given by Eqn. (B2). The IC photon num- 
ber flux Ty is derived by means of volume integration over the 
emission region and correct accounting for the growth of the area 
of the emission sphere on which the photons are distributed: 

l+z 



4nD^ 



dVAic[{l+z)Ei,{l+z)E2l 



(B13) 
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Here D denotes the luminosity distance and the additional factors 
of 1 + z account for the cosmological redshift of the photons. 



B4 7 -ray emission from decaying pions 

Provided the CR population has a power-law spectrum, the inte- 
grated y-ray source density Xy for pion decay induced y-rays can be 
obtained by integrating the y-ray source function Sy{Ey) (cf. EnBlin 
et al. 2007), 



Ay = Ay(El,E2) 



dEy Sy(Ey) 



4Cp m^oC(TppWN / rn^ V 



3a6y 



V 

2,26^ 



la+l a - I 

\ 26y ^ 26y 



1 + 



for /G{1,2}, 



(B14) 



,(B15) 



(B16) 



Cp is the normalisation 



where we used the abbreviation a = 
of the proton distribution function which we follow dynamically 
in our simulations (cf. Eqn. (A39)), and the rest mass of a neutral 
pion is m„oc^ ^135 MeV. The shape parameter 6y depends on the 
spectral index of the y-ray spectrum a according to 

6y^0.Ua-^-^ + 0A4. (B17) 

There is a detailed discussion in Pfrommer & EnBlin (2004a) how 
the y-ray spectral index ay relates to the spectral index of the parent 
CR population a^. In Dermer's model, the pion multiplicity is inde- 
pendent of energy yielding the relation ay = a (Dermer 1986a,b). 
The formalism underlying Eqns. (B14) and (B15) includes the de- 
tailed physical processes at the threshold of pion production like 
the velocity distribution of CRs, momentum dependent inelastic 
CR-proton cross section, and kaon decay channels. The y-ray lu- 
minosity is defined by 



£^ = JdVAy^Ay JdVCpH^, 



(B18) 



where we introduced the constant Ay with the dimensions [Ay] = 
y cm^ s~^ that is given by Ay = Ay/iC^n^) = const, according to 
Eqn. (B15). ^ is derived by Eqn. (B13) substituting Aic by Ay. 



It turns out that the primary synchrotron/IC emissivities scale as 
7v//c,prim oc leading to a scaling of the surface brightness of 
'S'v//c,prim ^ h^. In contrast, the secondary synchrotron/IC/y-ray 
emissivities scale as jV/zcsec ^ which results in a scaling of 
the surface brightness of 5v/7c,prim ^ h^- The diff'erent scaling of 
the primary and secondary non-thermal emission components with 
the Hubble constant is the reason why we choose to show all non- 
thermal luminosities in units of the currently favoured Hubble con- 
stant, hjo, where Hq = lOhjo km s~^ Mpc~^ 

This paper has been typeset from a T^/ MrgX file prepared by the 
author. 



B5 SPH projections and Hubble scaling 

We produced projected maps of the density, Mach number of 
shocks, relative CR pressure of protons and electrons, and non- 
thermal cluster observables in the radio, hard X-ray, and y-ray 
regime. Generally, a three-dimensional scalar field a(r) along any 
ray was calculated by distributing the product of a(r) and the spe- 
cific volume Mai pa of the gas particles over a grid comoving with 
the cosmic expansion. This yields the projected quantity A(r_L): 

A(r^,/;) = 7^ 2 ^c.,/;(r^,/; - rj, (B19) 

^pix a 

where Wajj is the value of the projected smoothing kernel (nor- 
malised to unity for the pixels covered) of an SPH particle a 
at comoving grid position r^jj, and Lp.^ is the comoving area 
of the pixel. In order to obtain a line-of-sight average of some 
mass-weighted quantity, say temperature, we project the quantity 
= TaPa dividcd each pixel by the mass projection (e.g. setting 

da = Pa)' 

Combining primary and secondary non-thermal emissivities 
requires the knowledge of the scaling with the Hubble constant. 
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